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1.  INTRODUCTION 


This  report  documents  tasks  performed  by  JAYCOR  from  the  period  June  1980 
until  October  1981  to  develop  a  biomechanical  understanding  of  the  processes 
Involved  in  air  blast  injury.  The  thrust  of  the  work  has  been  not  only  to 
carry  out  specific  tasks  but  to  evaluate  existing  technology,  test  and  digest 
alternate  approaches,  and  formulate  an  overall  plan  for  accomplishing  the 
mission  of  the  Blast  Overpressure  Program.  In  the  course  of  carrying  out  both 
the  spirit  and  the  letter  of  the  program  objectives,  we  have  investigated  a 
wide  range  of  phenomena  and  technology  that  had  been  previously  unknown  or  not 
applied  in  the  biomechanical  field.  The  results  have  been  encouraging  and  have 
allowed  the  Walter  Reed  Army  Institute  of  Research  (WRAIR)  to  be  able  to 
formulate  a  long-term  research  program  to  quantify  the  processes  and  develop 
new  damage  risk  criteria. 

The  project  during  this  period  of  time  consisted  of  three  contract  phases 
reflecting  the  changing  and  increasing  knowledge  about  the  mechanics  of  air 
blast  injury.  With  each  modification,  new  directions  that  had  been  explored  in 
the  previous  work  that  proved  promising  were  called  out  for  more  detailed 
investigation,  while  approaches  that  had  been  superseded  or  found  to  be  un¬ 
necessary  were  dropped.  Consequently,  varying  amount  of  work  was  done  on  each 
task  depending  on  the  viability  of  the  approach  taken.  In  doing  this  explora¬ 
tory  development,  close  contact  was  kept  with  WRAIR  and  the  decisions  to 
change  the  level  of  effort  on  the  major  tasks  was  done  with  their  knowledge 
and  direction.  The  result  has  been  a  highly  flexible  and  effective  working 
arrangement  between  sponsor  and  contractor  that  is  vital  to  a  scientific  in¬ 
vestigation  in  an  area  where  all  of  the  technical  issues  have  not  been  clearly 
identified. 

This  report  documents  many  of  the  details  of  the  work  during  the  period 
of  performance  grouped  by  their  scientifically  logical  order.  There  are  four 
main  parts,  Sections  2,  3,  4,  and  5,  that  address  in  causal  order  the  propa¬ 
gation  of  the  blast  wave,  the  interaction  with  a  body  shape,  the  state  of 
understanding  of  the  biomechanics  of  the  thorax,  and  finite  representations  of 
the  biomechanical  dynamics.  This  framework  contains  the  building  blocks  upon 


which  future  long-term  research  to  develop  and  validate  a  damage  risk  criteria 
for  air  blast  exposure  will  be  built. 


Contractually,  the  tasks  were  defined  during  the  course  of  the  work  and 
reflected  the  current  understanding  of  the  problem  and  therefore  included  work 
which  was  later  decided  to  be  expanded  upon  or  significantly  reduced.  The 
following  is  a  list  of  those  tasks  and  a  brief  summary  of  the  work  performed 
on  each.  There  were  a  total  of  eight  tasks  over  the  contract  period,  and  they 
have  been  renumbered  for  simplicity  but  appear  in  the  same  order  as  the  three 
contract  modifications. 

Task  1  -  Blast  Field  Interpolation 

The  BLAST  Code  was  used  to  Interpolate  to  all  points  in  the  field  the 
data  collected  from  M198  firings.  The  Intention  was  to  validate  the  calcula¬ 
tion  and  highlight  worst  case  pressures  distributions  around  the  weapon.  The 
comparisons  were  carried  out  in  detail  and  displayed  in  Section  2.1.  A  special 
contour  graphic  package  was  developed  for  displaying  the  output  of  the  com¬ 
puter  code  at  all  spaclal  locations  for  quantities  of  interest  such  as  peak 
overpressure,  A-duration,  and  A-impulse.  Static  as  well  as  dynamic  pressures 
were  determined  to  assess  the  Influence  of  the  winds  generated  by  the  explo¬ 
sions.  The  results  were  presented  to  WRAIR  in  an  Interim  progress  report. 

Task  2  ~  Similarity  of  Blast  From  Various  Sources 

The  purpose  of  this  task  was  to  evaluate  the  feasibility  of  using  the 
BLAST  Code  to  simulate  blasts  from  sources  other  than  the  gun.  In  particular, 
from  a  shock  tube  and  bare  charges.  At  the  time  the  task  was  written,  it  was 
uncertain  what  would  be  the  best  source  of  blast  waves  for  conducting  syste¬ 
matic  testing.  The  Impetus  for  this  comparison  was  diminished  when  the  pro¬ 
grammatic  decision  was  made  to  go  with  an  lmpactor  in  a  laboratory  environment 
as  providing  the  best  facilities  for  conducting  the  medical  research.  JAYCOR 
did,  however,  determine  that  the  BLAST  Code  would  be  applicable  to  bare 
charges  in  a  straightforward  modification  of  the  present  version  and  to  a 
shock  tube  when  calibrated  as  for  the  weapon  itself.  The  bare  charge  feature 
was  Implemented  as  an  option  to  the  code  and  a  description  of  its  use  and 
results  is  contained  in  Section  2.3. 


Task  3  -  Literature  Search  on  Biomechanical  Aspects 


JAYCOR  had  already  begun  a  literature  search  on  its  ovn  into  lung  and 
chest  models  and  this  task  extended  that  search  to  Include  data  bases  and  the 
thorax  region  as  yell  as  searching  for  lung  injury  mechanisms  directly.  The 
literature  searches  turned  up  almost  400  citations  which  were  relevant  to  the 
overall  problem,  and  these  citations  were  further  investigated  and  SO-some 
references  extremely  pertinent  to  the  mission  of  the  blast  overpressure  pro- 
gram  were  obtained  and  delivered  to  WRAIR.  In  addition,  the  complete  work  of 
Clemedson  on  the  subject  was  collected  and  delivered  in  a  bound  volume.  The 
results  of  that  search  are  described  in  Section  4.1,  and  Appendices  A  and  B 
contain  the  citations. 


Task  4  -  Conduct  a  Workshop  on  Biomechanical  Modeling 


Based  on  the  literature  search,  discussions  with  leading  experts  in  the 
biomechanical  modeling  field,  and  analysis  by  the  blast  overpressure  program 
manager,  it  was  decided  to  have  the  experts  meet  in  a  workshop  environment  in 
Albuquerque  during  December  1980.  JAYCOR  organized  and  hosted  the  workshop  and 
supplied  support  services  as  well  as  honoraria  to  the  guests.  Tape  recordings 
and  extensive  notes  were  taken  of  the  sessions,  and  after  the  workshop  a  syn¬ 
optic  report  on  the  presentations  and  conclusions  was  prepared.  That  material 
and  the  description  of  the  activities  is  contained  in  Section  4.2. 


Task  5  -  Develop  an  Three-Dimensional  Finite  Element  Representation  of 
Appropriate  Body  Shapes 


This  task  was  set  before  the  workshop  was  held  and  was  intended  to  begin 
the  process  of  building  a  finite  representation  of  the  body  for  both  external 
fluid  dynamics  calculations  and  internal  body  structural  calculations.  The 
task  would  draw  from  JAYCOR's  experience  in  both  finite  difference  and  finite 
element  modeling.  As  a  result  of  the  workshop,  however,  existing  modeling  of 
the  structural  aspects  were  revealed  and  it  was  the  consensus  of  both  WRAIR 
and  JAYCOR  that  two-dimensional  models  of  thorax  cross  section  should  be  in¬ 
vestigated  first.  JAYCOR  therefore  acquired  the  computer  programs  that  were 
relevant  to  this  task  and  began  adapting  them  to  the  blast  overpressure  mis¬ 
sion.  We  also  implemented  the  existing  lump  parameter  models  with  an  accurate 
solution  technique  into  a  convenient  computer  program  which  was  delivered  to 


WRAIR  for  its  in-house  use.  Furthermore,  in  discussions  with  Professor  Fung  at 
the  University  of  California,  San  Diego,  it  became  clear  that  lung  damage 
mechanisms  may  be  governed  by  wave  propagation  phenomena  within  the  lung 
parenchyma.  Therefore,  JAYCOR's  SPUNG  Code  was  adapted  to  a  body-lung  configu¬ 
ration  and  scoping  studies  made  on  the  feasibility  of  following  the  detailed 
wave  motion.  The  results  of  the  scoping  studies  are  contained  in  Section  5  and 
have  proven  to  be  important  in  focusing  the  effort  of  the  long-range  blast 
overpressure  program. 

Task  6  -  Further  Far  Field  Analysis 

As  a  result  of  the  interim  report  on  interpolating  the  M198  blast  field, 
WRAIR  personnel  identify  apparent  inconsistencies  in  some  far  field  values  of 
A-impulse.  The  calculations  were  reviewed  and  it  was  discovered  that  calibra¬ 
tion  based  on  Initial  peak  overpressure  was  subject  to  considerable  error  be¬ 
cause  of  the  noisiness  of  the  signal.  In  several  cases,  we  had  chosen  an 
overly  large  value  at  nearby  points  for  calibration  which  invalidated  the  far 
field  results.  At  WRAIR's  suggestion,  new  calculations  were  made  based  on  cal¬ 
ibrating  on  A-impulse,  an  Integral  quantity  that  is  expected  to  be  less  sensi¬ 
tive  to  the  signal  noise.  That  method  of  calibration  has  indeed  proved  to  be 
effective  and  the  comparisons  are  shown  in  Section  2.1.  We  now  are  able  to 
reproduce  not  only  the  qualitative  nature  of  the  pressure  signal  but  the  quan¬ 
titative  variation  with  distance  from  the  source  of  A-duration,  A-impulse,  and 
peak  maximum  pressure. 

Task  7  -  Near  Field  Blast  Modeling 

The  purpose  of  this  task  was  to  make  detailed  calculations  of  the  loading 
on  a  two-dimensional  object  placed  in  a  blast  field  and  to  prepare  a  protocol 
for  carrying  out  testing  that  would  validate  the  results.  Three  shapes  were 
identified,  ellipse,  square,  and  circle.  Calculations  were  made  for  each  case 
with  a  3  psi  and  a  20  psi  blast  wave.  Two  orientations  of  the  ellipse  were 
investigated  also.  The  results  are  discussed  in  Section  3.1. 


Because  of  higher  priority  testing  programs,  Lovelace  Research  Institute 
which  was  to  carry  out  the  testing  phase  has  been  unable  to  make  testing  time 


available  during  the  current  contract  period.  Therefore,  under  WRAlR's  direc¬ 
tion,  we  formulated  a  protocol  and  preliminary  design  of  a  test  target.  The 
protocol  is  found  In  Section  3.2,  but  the  experiment  will  have  to  be  carried 
out  later  in  FY  1982  at  WRAlR's  discretion. 

Task  8  -  Modeline  of  Bodv  and  Lune  Dynamics 


Because  of  the  encouraging  results  shown  by  General  Motors  at  the  bio¬ 
mechanical  workshop  and  those  developed  by  JAYCOR  under  Task  5,  it  was  decided 
to  intensify  the  effort  to  use  the  existing  structural  analysis  program  called 
FEAP  (Finite  Element  Analysis  Program).  At  first  it  was  believed  that  the  code 
could  be  acquired  and  used  immediately,  but  because  of  the  proprietary  nature 
of  General  Motors'  research,  JAYCOR  was  forced  to  work  with  the  developers  of 
the  code  at  the  "niversity  of  California,  Berkeley,  to  reconstruct  the  analy¬ 
sis.  In  doing  so,  we  uncovered  and  solved  coding  problems  in  the  program  and 
made  corrections  to  the  material  properties  that  had  been  suggested  by  General 
Motors.  The  effort  to  prepare  the  code  for  use  in  the  blast  overpressure  pro¬ 
gram  required  more  effort  than  had  originally  been  expected.  However,  we  have 
now  overcome  those  difficulties  and  can  use  the  program  in  both  a  static  and 
dynamic  mode.  Section  5.1  describes  the  first  analyses  using  FEAP  and  this 
effort  is  continuing  in  anticipation  of  its  application  in  the  next  phase  of 
the  blast  overpressure  program. 


In  summary,  the  tasks  performed  under  this  contract  and  its  modifications 
have  scoped  and  developed  the  technologies  required  to  assist  the  blast  over¬ 
pressure  program  in  quantifying  and  validating  the  damage  risk  criteria  for 
future  weapon  deployment.  As  was  mentioned  earlier,  this  effort  has  been 
directed  toward  guiding  the  use  of  technology  as  well  as  implementing  it  into 
a  specific  form.  The  availability  of  complete  and  consistent  models  connecting 
the  blast  source  through  the  far  field  propagation  to  the  loading  on  the  body 
to  the  structural  response,  and  finally  to  the  local  stress  distributions 
within  the  lung  is  an  important  part  of  the  blast  overpressure  program 
mission. 


2.  BLAST  OVERPRESSURE  FAR  FIELD 


2.1  THEORY  OF  BLAST  CODE 


2.1.1  Blast  Have 


A  blast  overpressure  wave  Is  created  when  high  pressure  and  temperature 
product  gases  that  propel  the  shell  leave  the  gun  barrel.  The  disturbance  to 
the  atmosphere  changes  rapidly  in  amplitude  and  shape  as  it  propagates*  At  the 
front  of  the  blast  wave,  the  pressure  p  and  density  p  jump  abruptly  from  their 
undisturbed  values.  Immediately  after  the  front  passes,  this  disturbance  re¬ 
turns  to  ambient  quickly  and  is  followed  by  a  rarefaction  wave  (Figure  1).  The 
structure  of  the  blast  wave,  e.g.,  the  amplitude  and  the  duration  of  both  the 
compression  and  the  rarefaction  parts,  varies  according  to  its  source  strength 
and  ambient  conditions.  In  the  standard  atmosphere,  an  extremely  strong  source 
such  as  a  nuclear  explosion  will  generate  a  strong  blast  wave,  characterized 
by  a  shock  front  and  a  very  shallow  rarefaction  tail.  On  the  other  hand,  a 
weak  source  will  generate  a  sound  wave  whose  amplitude  and  the  duration  are 
equal  for  the  compression  and  the  rarefaction  parts  of  the  wave.  A  blast  wave 
with  source  strength  between  these  two  extreme  cases  will  have  a  wave  form  of 
mixed  type.  The  variation  of  the  wave  form  with  the  source  strength  is  indi¬ 
cated  schematically  in  Figure  2. 

The  most  appropriate  mathematical  treatment  of  the  waves  also  varies  with 
the  strength  of  the  source.  Since  sound  waves  are  weak  disturbances,  the  per¬ 
turbation  to  the  flow  variables,  p',  p“  and  u',  are  small  quantities,  with  the 
perturbation  to  entropy,  s'*,  being  three  orders  of  magnitude  smaller  [Ref.  1]. 
Therefore,  the  entropy  of  the  sound  waves  is  essentially  constant.  The  gov¬ 
erning  equations  can  be  linearized  by  dropping  higher  order  terms  and  reduce 
to  a  simple  wave  equation.  The  technique  for  solving  this  equation  can  be 
found  in  standard  textbooks  (such  as  Ref.  2). 

To  study  the  opposite  extreme  of  large  source  strength,  the  concept  of 
similarity  was  introduced  independently  by  Sedov  [1946,  Ref.  3]  and  Taylor 
[1950,  Ref.  4],  The  concept  has  been  used  in  other  branches  of  fluid  dynamics, 


N 
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W.D.  =  Weak  discontinuity 
S.W.  *  Shock  wave 


Figure  2.  Pressure  Trace  Variation  with  Source  Strength. 
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such  as  boundary  layer  theory,  conical  flow  theory,  and  transonic  and  hyper¬ 
sonic  flow  theory.  The  assumption  of  similarity  decreases  the  number  of  inde¬ 
pendent  variables  and  often  reduces  the  governing  partial  differential 
equations  to  more  manageable  ordinary  differential  equations.  The  technique 
has  been  very  successful  in  analyzing  strong  blasts. 

For  blast  waves  with  intermediate  source  strength,  numerical  methods  have 
to  be  employed  to  solve  the  hyperbolic  gasdynamics  equations  and  one  of  the 
best  is  the  method  of  characteristics.  The  distinguishing  property  of  the 
hyperbolic  equations  is  the  existence  of  certain  special  directions  or  lines 
in  the  r-t  plane  called  characteristics  (see  Ref.  5).  Along  the  characteris¬ 
tics,  the  dependent  variables  satisfy  certain  equations  known  as  compatibility 
relations.  Solution  of  the  compatibility  relations  generates  the  space-time 
history  of  blast  wave,  however  the  procedure  is  complicated  since  the  charac¬ 
teristic  lines  are  themselves  unknowns  to  be  determined. 

The  overpressure  data  collected  for  the  M198  howitzer  with  a  M-203  charge 

[Ref.  6]  indicates  that  these  waves  are  of  intermediate  source  strength.  It 

can  be  shown  that  for  the  strong  blast  wave  the  attenuation  of  the  maximum 

pressure  is  r  while  for  the  sound  wave  it  is  r  [Ref.  2],  A  systematic 

analysis  (Table  1)  of  the  experimental  gun  overpressure  data  indicates  that 

—1  —2 

the  maximum  overpressure  attenuates  as  r  to  r  ,  indicating  that  the  waves 

are  of  intermediate  strength  and  requiring  the  use  of  numerical  techniques. 

For  instance,  Table  1  shows  the  measured  values  of  the  amplitude  of  the  inci- 

—1  —2 

dent  wave  and  its  comparison  with  the  values  obtained  from  the  ~r  and  ~r 
relationships.  The  BLAST  code  which  employs  the  method  of  characteristics  is 
developed  to  handle  blast  overpressure  waves  typical  of  those  generated  by  the 
M198  howitzer  and  similar  weapons. 

Table  l. 


Distance 
Defined  in 
Figure  1 

Distance  from 
the  Muzzle 
Brake,  r 

~  r'1 

Experimental 

(m) 

(m) 

(psig) 

(psig) 

10 

11.44 

1.41 

2.50 

,  r.  «r.  it;  v;  A  ■.  P,1  -■ 


■  »  ■  f "  >  ->  ■  -  '1  'I  'T 


2.1.2  Quasispherlcal  Approximation 

After  the  explosion  of  the  charge,  the  shell  is  propelled  by  the  high 
pressure  and  high  temperature  product  gas,  which  leaves  the  muzzle  brake  as 
the  shell  is  launched.  Initially,  the  gas  flow  in  a  small  region  surrounding 
the  muzzle  brake  is  complicated  by  the  geometries  of  the  gun  and  the  muzzle 
brake.  This  effect,  however,  becomes  less  significant  at  distances  large  com¬ 
pared  to  the  source  diameter  and  the  wave  front  propagates  more  like  a  spheri¬ 
cal  wave. 

In  order  to  take  advantage  of  the  simple,  near  spherical  structure  of  the 
blast  wave  at  large  distance,  the  BLAST  code  employs  a  "quasi-spherical"  de¬ 
scription.  That  is,  along  each  radial  line  emerging  from  the  end  of  the  gun 
barrel  the  wave  is  treated  as  though  it  is  part  of  a  purely  spherical  wave. 
The  intensity  of  the  equivalent  spherical  wave  is  allowed  to  vary  with  the 
direction,  so  that  the  effect  of  the  nonspherlcal  geometry  of  the  muzzle  brake 
can  be  simulated. 

In  this  quasi-spherical  approximation,  the  calculation  along  each  radial 
line  is  independent  from  all  other  radial  lines.  Along  each  radial  line  only  a 
spatial  coordinate,  i.e.,  the  radial  distance  from  the  source,  and  the  time 
variable  are  involved.  Thus  the  approximation  essentially  reduces  a  compli¬ 
cated  three  dimensional  problem  into  a  set  of  one  dimensional  problems.  The 
method  of  characteristics  for  one  dimensional  gas  dynamics  is  readily  applica¬ 
ble. 

Ground  reflection  is  an  important  part  of  the  blast  waves  in  M198  and 
similar  gun  firings.  Unfortunately  the  reflection  of  blast  wave  is  difficult 
to  handle  both  theoretically  and  numerically.  A  complete  treatment  of  the  re¬ 
flection  should  consider  the  ground  absorption  and  the  nonlinear  effects  such 
as  the  variation  of  the  reflection  angle  and  the  formation  of  the  Mach  stem. 
Such  a  treatment  is  far  too  complicated  to  be  investigated  in  the  scope  of  the 
present  project. 

Instead,  we  employ  a  simpler  approach  to  the  reflection  problem,  namely, 
we  neglect  the  ground  absorption  and  nonlinear  effects  and  treat  the  reflec¬ 
tion  as  perfect.  The  viability  of  this  approach  is  borne  out  in  the  comparison 
with  field  data. 
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For  perfect  reflection,  the  method  of  image  can  be  used.  Thus  the  re¬ 
flected  wave  is  regarded  as  generated  from  an  image  source  at  a  distance  below 
the  ground  level  equal  to  the  height  H  of  the  muzzle  brake  (see  Figure  3).  The 
pressure  wave  impinging  on  the  pressure  sensor  PS  is  then  a  linear  superposi¬ 
tion  of  the  direct  incident  wave  and  the  reflected  wave.  The  quasi-spherical 
approximation  and  the  method  of  characteristics  are  used  both  for  the  Incident 
and  reflected  waves. 

2.1.3  Pronaeation  of  Blast  Wave  Alone  a  Radial  Line 


We  choose  the  center  of  the  muzzle  brake  as  the  origin  of  a  spherical 
coordinate  system  (r,  0,  t)  shown  in  Figure  4.  Assuming  the  gas  is  ideal, 
inviscid  and  flowing  isentroplcally, *  the  equations  of  continuity,  motion,  and 
energy  can  be  written: 

_£.+  i_  1  +  .  9-  —  (pv  sin  0)  +  — sin_0  3*  (wp)  “  0  (1) 

r 

_|du  v2  +  w2l  _  3p 


Dv  .  uv  w  cot  0  1  3p 


Dw  uw  vw  cot  0 

'  Dt  r  r 


1  3£_ 

r  sin  0  3$ 


D  .  ,  2  D  .  . 

Dt  [p]  ’  a  Dt  [Pl 


where 


D  _  3__  v  _3_  w  3 

Dt  “  3t  u  3r  r  39  r  sin  6  3^ 


*The  isentropic  assumption  would  not  be  valid  for  shock  waves  associated  with 
8 1 rong  blasts  and  the  Rankine-Hugoniot  shock  relations  would  have  to  be 
employed  to  account  for  the  entropy  Jump.  The  assumption  is  valid  in  the 
present  case. 


Figure  3.  Method  of  Image 


Figure  A.  Spherical  Coordinates  with  Origin  at 
Center  of  Muzzle  Brake. 
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Here  u,  v  and  w  are  the  velocity  components,  p  Is  the  pressure,  p  is  the  den¬ 
sity,  a  is  the  local  sonic  speed,  and  t  is  the  time  measured  from  the  instant 
that  the  blast  waves  are  generated. 

A  complete  solution  to  Equations  (1)  to  (5)  depends  on  the  source  distri¬ 
bution  which  in  turn  depends  on  the  detonation  of  the  explosive  charge  and  the 
resulting  flow  field  in  the  barrel  and  nuzzle  brake.  In  this  work  we  shall 
assume  that  the  source  distribution  is  known  and  not  attempt  to  calculate  it 
from  more  fundamental  processes.  In  the  quasi-spherical  approximation  employed 
in  this  work,  we  simulate  the  origin  of  the  blast  waves  along  each  radial 
direction  as  a  sphere  of  pressurized  gas  and  the  transport  processes  are  im¬ 
portant  only  along  the  radial  direction  although  the  source  strength  may  vary 
slowly  with  6  and  i|>.  In  this  approximation.  Equations  (1)  to  (5)  can  be  sim¬ 
plified  as  follows: 

A 

Pfc  +  uPr  +  p  ur  +  -  0  ,  (6) 

ut  +  uur  +  £  pr  -  0  ,  (7) 

Pt  +  upr  “  a2(pt  +  upr ^  “  0  ♦  (8> 

where  the  subscript  t  and  r  denote  partial  differentiation. 

The  initial  and  boundary  conditions  are  taken  as 

At  t  «  O', 


for  r  <  r. 


P  -  P„ 


p  -  p 


u  ■  0 


for  r  >  r„ 


•  ^  ^  v*  cl  cf 4.^0 


where  rQ  *  initial  radius  of  the  pressurized  gas  sphere  and  the  subscripts  s 
and  o  denote  source  and  undisturbed  conditions,  respectively. 

At  t  >  0+, 


u  -  0 


at  r  -  0 


Since  Equations  (6)  to  (8)  are  hyperbolic  equations  [Ref.  8],  their  char¬ 
acteristic  form  can  be  obtained  in  the  following  way.  First,  let  us  define 

o  z  f  &  .  (12) 

J  pa 

Substitution  of  Equation  (12)  into  Equations  (6)  to  (7)  yields 


0  +  uo  +  au  +  -  0 

t  r  r  r 


u  +  uu  +  ao  ■  0  . 

t  r  r 


Adding  and  subtracting  Equations  (13)  and  (14),  we  obtain 


(a  +  u)t  +  (u  +  a)  (o  +  u)f  +  ■  0  , 

(o  -  u)t  ♦  (u  -  a)  (o  -  u)r  +  -  0  , 

respectively.  It  may  be  shorn  that  for  lsentroplc  flow 


0  -  — f  • 

Equations  (15),  (16)  and  (8),  with  the  substitution  of  Eq.  (17),  can  then 
be  cast  into  characteristic  form  as  follows 
On  f*’  curve: 


dr 

dt"  "  u  * 


d  /  2a  \  2au  n 

3f  [r-—  ° 


On  r"  curve: 


d  /  2a  \  .  2au  _  n 

ir  {rrr  "r  —  ■  0 
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(22) 


•.rv  c  %r 


On  r° 


iL  _  J 


where  d/d T  is  the  total  differentiation  with  respect  to  t  along  the  corre¬ 
sponding  characteristics. 

Equations  (18),  (20)  and  (22)  define  the  direction  of  the  characteris¬ 
tics  r+,  T  and  r°  respectively.  Along  these  characteristics,  the  compatibil¬ 
ity  relations,  i.e.,  Equations  (19),  (21)  and  (23),  are  satisfied  accordingly. 
A  sketch  of  T+,  T~  and  r®  characteristics  are  shown  in  Figure  5.  It  should  be 
noted  that  the  characteristics  r®  are  identical  to  the  streak  lines  of  the 
fluid  particle.  Furthermore,  the  blast  wave  front  follows  closely  with  one  of 
the  r+  characteristics. 

Separate  treatments  for  the  blast  wave  front  and  the  other  part  of  the 
blast  wave  are  necessary.  At  the  wave  front  the  pressure,  density  and  particle 
velocity  are  discontinuous.  Such  discontinuities  are  described  by  the  shock 
adiabatics  [Ref.  7] 


1  a2  (l±±  P_  +  T  -  l\ 

2  a°  \  y  PG  V  )  ' 

4(H  • 


(Y  +  l)p  +  (Y  -  l)p0 

P  "  Po  (Y  -  l)p  +  (Y  +  1)P0  ’  (26) 

where  ufl^  is  the  velocity  of  the  shock  front  and  p,  u  and  p  are  the  pressure, 
fluid  velocity  and  density  just  behind  the  shock  front,  respectively.  Notice 
that  through  Equations  (24)-(26)  all  the  four  quantities  ugh,  p,  u  and  p  are 
determined  when  one  of  them,  ug^  say,  is  known.  When  p  and  p  are  known,  the 
velocity  of  sound  just  behind  the  shock  front  can  also  be  calculated  from  the 


SS'j'  .v.NV. 
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Figure  5.  Characteristics  In  r-t  Plane. 

relation  a?"  m  Yp/p.  Now  the  pair  of  equations  (18)  and  (19)  of  the  T+  charac¬ 
teristics  reduces  essentially  to  involve  only  one  unknown  flow  variable  and 
therefore  can  be  solved. 

For  other  parts  of  the  wave,  all  three  characteristics  r+,  r~  and  T®  are 
needed  to  determine  the  flow  variables.  Equations  (18)  to  (21)  consist  of  only 
two  unknown  flow  variables,  u  and  a.  These  equations  are  solved  first.  Once  u 
and  a  are  found.  Equations  (22)  and  (23)  can  be  used  to  obtain  the  solution  of 
p  after  the  substitution  of  p  ■  Yp/a  . 

The  solution  of  the  characteristic  equations,  Eqs.  (18)  to  (21),  is  based 
on  the  upstream  interpolation  scheme  introduced  by  Belotserkovskii  and  Chush- 
kin  [Ref.  8].  The  basic  idea  is  as  follows:  Firstly,  approximate  those  equa¬ 
tions  by  use  of  first  order  implicit  finite  difference  formula  for  small  time 
steps,  which  results  in  a  piecewise  linear  characteristic  network.  Implicit 
finite  difference  method  has  the  advantage  that  it  is  stable  and  therefore 
allows  the  use  of  larger  sizes  of  dr  and  At.  The  wave  front  is  traced  by  fol¬ 
lowing  one  of  the  T+  characteristics,  initiated  at  the  boundary  of  the  pres¬ 
surized  gas  sphere.  This  particular  characteristic  divides  the  disturbed  flow 
region  from  the  undisturbed  one.  A  number  of  points  in  r-direction  with  equal 
spacing  6r  are  used.  The  exact  number  depends  on  how  far  from  the  muzzle  brake 
we  want  to  calculate.  Convergent  solutions  at  one  time  level  are  extrapolated 
to  give  initial  guesses  for  flow  variables  at  the  next  time  level.  They  are 
substituted  into  Equations  (18)  and  (20)  which  give  roughly  the  locations  of 
the  r+  and  T“  curves  at  old  time  level.  Hence  flow  variables  at  these  loca¬ 
tions  can  be  determined  by  Interpolation  between  the  convergent  solutions  at 
the  old  time  level.  With  this  information  fresh  values  at  new  time  level  are 
obtained  from  Equations  (19)  and  (21).  These  fresh  values  at  new  time  level 


can  be  substituted  into  Equations  (18)  and  (20)  again  to  initiate  another 
cycle  of  iteration.  This  iteration  process  is  repeated  until  two  successive 
guesses  of  the  flow  variables  at  the  new  time  level  agree  within  sufficiently 
close  limits.  Due  to  the  implicit  nature  of  the  method,  similar  procedure  is 
used  in  calculating  the  wave  front  and  in  solving  p  and  p  in  Equations  (22) 
and  (23).  The  solution  procedure  can  be  summarized  in  the  following  flow  chart 
(Figure  6).  More  detailed  discussion  of  the  iteration  procedure  is  presented 
in  Section  2.1.4. 

2.1.4  Iteration  Procedure  for  the  r+  and  T~  Characteristics 

We  present  the  iteration  procedure  for  the  and  T  characteristics  in 
detail  in  the  following.  The  iteration  procedures  on  the  T®  characteristics 
and  for  the  shock  front  are  similar  and  are  not  detailed  here. 

Let  r^,  i  ■  1,  ...,  N  be  a  set  of  equally  spaced  mesh  points  of  the 
radial  variable  r,  which  is  fixed  once  of  all,  for  all  time  steps.  We  consider 
the  finite  difference  equations  of  (18)-(21)  between  the  instants  t  and  t  +  fit 
in  the  following  form: 

ri  “  rI±(t)  *  t]  +  u[r^ ,  t  +  fit]  ±  a[r1±(t),  t]  +  a^,  t  +  fit]} 

(27) 

J±<ri»  t  +  6t>  "  J±lri±(t),  tJ  "  “  f^-  Mr1±(t),  t]  +  f[rt,  t  +  fit  ] }  (28) 

where 

^(r.t)  -  ±  u(r,t)  (29) 

f(r,t)  -  2a(r,t)  u(r,t)/r  (30) 

In  these  equations,  the  alternative  signs  ±  are  for  T+  or  r“  characteristics, 
respectively.  The  quantities  r^(t)  are  the  points  on  the  T*  characteristics, 
which  will  propagate  to  r^  at  t  +  fit.  Note  that  r^(t)  do  not  in  general 
locate  at  the  mesh  points  r^.  We  assume  that  the  flow  variables  at  the  mesh 


**The  superscript  v  denotes  the  vth  Iteration,  and  e  is  the  tolerance  limit  of 
relative  error  for  convergence. 


Figure  6.  Flow  Chart  for  the  Method  of  Characteristics. 


point  at  time  t  [i.e.,  u(r,t)  and  a(r,t)]  and  those  at  previous  time  steps 
are  known.  We  need  to  find  the  flow  variables  at  the  mesh  points  at  time 
t  +  fit  [i.e.,  u(r^,t+5t)  and  a(r£,t+fit)]. 

Equations  (27)-(28)  are  implicit  finite  difference  equations  because  the 
quantities  multiplied  by  fit  in  the  right  hand  side  of  these  equations  involve 
the  unknown  flow  variables.  These  equations  also  involve  the  unknown  r^±(t). 
We  solve  these  equations  by  the  following  Iteration  procedure. 

Let  us  denote  the  initial  guess  with  a  superscript  0  and  the  subsequent 

Iterated  values  with  superscripts  v  *  1,  2 . The  initial  guess  u®(r^, 

t+6t)  is  obtained  for  the  fluid  velocity  from  an  extrapolation  using  the 
values  u(r£,  t)  and  u(r^,  t  -  fit).  A  similar  initial  guess  a®(r^,  t+5t)  is 
obtained  for  the  velocity  of  sound.  The  initial  values  r®±(t)  of  r1±(t)  are 
then  calculated  from 

r°±(t)  -  r±  -  {u°(r1,  t  +  fit)  ±  a°(r1,  t  +  fit)}  .  (31) 


With  this  set  of  initial  values,  we  can  start  the  iteration. 

In  the  vth  cycle  of  the  iteration  uv-1(ri,t+fit),  av-1(r1,t+5t)  and 
r^^(t)  are  known  and  we  need  to  determine  their  new  iteration  values 
uv(r1,t+5t),  av(r^,t+fit)  and  r}^(t).  This  is  executed  as  follows.  First, 
u[r^1(t),t]  and  a[r]^(t),tj  are  calculated  by  interpolation  from  the  set  of 
u( r j , t )  and  a(rj,t),  j  ■  1,  ...,  N.  We  can  then  calculate  r)*±(t)  from 


Now  using  Equation  (29),  we  have 


uV(r4,  t+6t)  -  5-  {J^(r± ,  t+St)  +  t+6t)}  (34) 

aV(>r1»  t+it)  -  —  j  1  {J^(r1,  t+6t)  -  J^r^  t+«t)}  (35) 

and  these  are  the  new  Iteration  values  we  sought. 

The  iteration  is  repeated  until  the  following  convergence  conditions  are 
fulfilled: 

|{«V(r1»  t+6t)  -  uV  1(ri,  t+St}}/uV|  <  cy  (36) 

|{aV(r1,  t+6t)  -  aV“1(r1,  t+6t)}/aV|  <  c&  (37) 

where  eu  and  efl  are  predetermined  small  quantities.  When  these  conditions  are 
satisfied,  one  can  see  that  the  Equations  (32)  and  (33)  are  approximations  of 
the  original  finite  difference  equations  (27)  and  (28),  with  rti(t)  ■  ri±^» 
u(rif  t+6t)  -  uv(rit  t+$t)  and  a(r^,  t+6t)  ■  av(ri,  t+6t)  as  the  approximate 
solutions. 

2.1.5  Input -Output  of  the  BLAST  Code 

The  detailed  input  to  the  BLAST  code  is  described  in  Table  2.  However, 
most  of  the  parameters  there  are  related  to  the  numerical  methods  and  output 
options.  Therefore  Table  3  Indicates  a  shorter  list  that  can  be  used  in 
conjunction  with  the  default  values. 

The  code  presents  pressure-time  histories  and  other  Important  physical 
quantities  at  specified  spatial  locations,  in  the  form  of  tables  or  graphics. 
Some  of  these  quantities  are  as  follows: 

(1)  P^:  the  direct  incident  overpressure  in  psi. 

(2)  ?£:  the  reflected  component  of  the  overpressure  wave,  in  psi. 

Note  that  the  maximum  of  P2  is  not  equal  to  the  height  of  the 
second  peak  in  the  pressure-time  history. 


Card 

FORTRAN 

Identifier 

Meaning 

Units 

CARD  1 

DR 

Spatial  Increment  In  radial  direction 

ft 

DT 

Time  step  of  calculation 

sec 

TMAX 

Maximum  time  of  problem 

sec 

IMAX 

Maximum  number  of  spatial  points  along 

— 

radial  line 

CARD  2 

RF 

Radius  of  spherical  source 

ft 

PO 

Ambient  atmospheric  pressure 

psi  , 

RHO 

Ambient  atmospheric  density 

8lug/ftJ 

GAMMA 

Ratio  of  specific  heats 

— — 

CARD  3 

11 

Number  of  radial  points  to  be  contoured 

— 

KK 

Number  of  different  gun  elevations 

CARD  4 

HPIVOT 

Height  of  gun  pivot  point 

ft 

(total 

GUNLENF 

Length  of  gun  barrel  from  pivot  point 

ft 

of  KK 

ELEVATION 

Elevation  angle 

deg 

cards ) 

CARD  5 

IFLALNG 

Flag  to  obtain  lung  calculation 

— 

IFLGP 

Flag  to  obtain  plot  graphics 

— 

IFLGTAP 

Flag  to  obtain  tape  storage 

— 

CARD  6 

NTOT 

Total  number  of  radial  lines 

— 

NP 

Number  of  radial  lines  to  be  contoured 

— 

CARD  7 

PSI(I), 

Angular  position  of  rays 

deg 

PS  I (NTOT) 

CARD  8 

II 

Total  number  of  radial  points  to  be 

— 

calculated. 

CARD  9 

X(I1+1), 

Coordinates  of  points  to  be  calculated  but 

ft 

...»  X(II) 

not  plotted. 

CARD  10 

PS 

Source  pressure 

atm 

One  set 

of  cards  #8,  9, 

10  for  each  radial  line. 
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(3) 

Pstat: 

the  peak  static  overpressure,  in  psi.  It 

of  the  overpressure  in  the  pressure-time 

is  the  maximum  value 

history. 

(4) 

Pdyn: 

peak  dynamic  pressure,  in  psi. 

(5) 

A  • 

Adur* 

A-duration,  in  seconds. 

(6) 

^imp: 

A-impulse,  in  psi-sec. 

Table  3.  Short  Form  of  Input. 

I.  Ocher  elevations  of  previously  calibrated  gun. 
ELEVATN  Gun  Elevation 

II.  Different  muzzle  brake  design. 

PS(1) . PS(NP)  Source  strength  along  each  ray 
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As  a  validation  of  the  BLAST  code,  we  compare  here  results  of  the  BLAST 
code  calculation  with  the  M198  howitzer  data  in  the  May  15,  1979  firings,  as 
tabulated  in  the  Final  Report  prepared  by  JAYCOR  under  the  contract  with  WRAIR 
(No.  DAMD17-78-C-8087). 

Figure  7  shows  the  ground  map  for  all  the  locations  where  pressure  mea¬ 
surements  were  taken  and  at  which  comparison  with  calculation  is  made.  In  the 
measurement,  the  gun  was  set  at  three  elevation  angles,  equal  to  800,  267  and 
45  mils  (45°,  15°  and  2.53°,  respectively)  and  was  fired  several  times  so  that 
for  each  elevation  angle  and  at  each  location  three  sets  of  data  were  record¬ 
ed.  Part  of  the  data  is  reproduced  together  with  results  of  our  calculation  in 
Figures  9-11  and  16-18. 

As  a  first  test  of  the  BLAST  code,  calculations  are  performed  by  choosing 
a  set  of  initial  pressures  P8,  one  for  each  radial  line  (at  azimuthal  direc¬ 
tions  ■  0°,  30®,  60°,  120®  and  150°),  so  that  the  calculated  and  the  mea¬ 
sured  maximum  overpressure  match  each  other  at  a  near  field  location  (10 
meter).  The  Pfl  values  are  shown  In  Figure  8  for  the  three  elevation  angles 
considered.  The  initial  pressurized  balloon  radius  rQ  used  in  this  calculation 
is  3  feet. 

The  results  of  the  calculation  have  been  reported  in  the  JAYCOR  progress 
report  "Far-Field  Model  and  Validation,"  and  they  are  reproduced  here  in  Fig¬ 
ures  9-14.  In  Figures  9-11  the  calculated  and  the  measured  pressure-time  his¬ 
tories  are  compared  side  by  side,  and  In  Figures  12-14  the  distribution  of 
various  quantities,  such  as  the  maximum  overpressure,  A-impulse  and 
A-duration,  on  all  field  locations  is  presented  In  the  form  of  contour  graphs. 

In  Figures  9-11  one  can  see  that  the  calculation  has  reproduced  the  main 
features  of  the  measured  pressure  trace,  such  as  the  sharp  rise  of  the  pres¬ 
sure  at  the  shock  front,  the  arrival  times  of  the  direct  incident  and  ground 
reflection  peaks,  the  rapid  transition  from  compression  and  rarefaction  and 
the  final,  slow  relaxation  of  the  rarefaction  to  ambient.  This  shows  that  the 
method  of  characteristics  is  basically  a  correct  approach  to  the  problem. 


igure  8.  Source  Pressures  as  Calibrated  with  Peak  Pressure  for  the  May  1979  Firings 


mo 


MIC  I  P  PSIC  I  I  p  PSIC 


9.0  HH 


.95  9.19  9.  IS 

TIIC  SEC 
P  TRACE 


9.99  0.95  9.10  0.15  0.29 

■  IME  SEC 

p  TRACE 


- 

•  | 

) 

i 

^  ! 

a. 

0.0  - - '• 

i 

\  •- 

V 

3.00 

0.05  0.10  0.16  3.20 

TIME  SEC 

P  TRACE 

9.05  9.10  0.15 

T I  HE  SEC 
P  TRACE 


•'  •  «V  0.0*  0.10  9.15 

• ; Mr  EC 

P  TRACE 


>^\a>:vv 


PSIC 


PSIG 


0.  10 

Tl«  SEC 
P  TRACE 


0.  10 
TIME  SEC 
P  TRACE 


2.0  • 

(  i 

i 

o 

: 

£ 

\ 

a  a  -  ■-  \ 

i 

0. 

0.00  0.05 

0.10  0.IS 

0.20 

TIME  SEC 

P  TRACE 

0.00 

0.05 

0. 10 

0.  IS 

0.20 

TIME  SEC 

P  TRACE 

TIME  SEC 
P  TRACE 


PSIG  I  I  P  PSIG  !  I  P  PSIG  I  I  P  PS1G 


&  r/' 


PSIG  I  I  P  PSIG  I  I  P  PSIG 


2.0 


0.00 

0.05 

0.10 

0.15 

0.20 

TIME  SEC 

P  TRACE 

0.00 

0.05 

0.  10 

0.15 

0.20 

TIME  SEC 

P  TRACE 

0. 10 

tiw  sec 

P  TRACE 


0.  10 
TIME  SEC 
P  TRACE 


TIME  SEC 
P  TRACE 


TIME  SEC 
P  TRACE 


0.  10 
TIME  SEC 
P  TRACE 


PSIC  I  I  P  PSIC  II  p  PS,G 


SISd  4  I  I  SIM  *  I  I  SIM  4  I  I  SIM 


■  -  • 

-  •  ■  •  . - 

|| 

•3 

sn 

to 

1  :• : 

CL 

i !' 

0,0  - 1  1  ^  _ 

0. 

_o  a  j  _____  i _ 

c .*  _____  * 

0.00  0.05 

0. 10  0.15  0.20 

T IIC  SEC 

P  TRACE 

0.10 

TIME  SEC 
t  TRACE 


PSIG 


&  . 

nimp 


» 


V;V  \  VN 


vs  i 


A  closer  examination  reveals  that  certain  peaks  In  the  calculated  trace 
(e.g.,  the  traces  at  30  and  40  meters  In  Figures  9a-g)  are  too  large  compared 
with  the  corresponding  measured  values.  However,  one  can  also  see  in  these 
figures  that  there  are  sometimes  sharp  spikes  superimposing  on  the  measured 
trace,  especially  at  the  shock  fronts  at  the  near  locations  (10  and  20 
meters).  The  cause  of  these  spikes  is  not  clear  at  the  moment,  but  they  affect 
strongly  the  maximum  pressure  which  we  used  to  determine  the  initial  condi¬ 
tion,  i.e.,  the  initial  pressure  Pg  in  our  calculation.  It  is  preferable  to 
build  our  calculation  on  a  physical  quantity  which  does  not  depend  so  much  on 
these  spikes  of  unknown  cause  and  then  to  make  detail  quantitative  comparison 
with  the  field  data.  Such  a  physical  quantity  is  the  A-impulse,  since  it  is  an 
integrated  quantity  and  is  insensitive  to  rapid  fluctuations  and  narrow 
spikes.  We  present  the  results  of  such  a  calculation  in  the  next  subsection 
(Section  2.2.2). 

2.2.2  Comparison  by  Matching  A-Impulse 

For  a  fixed  elevation  we  choose  for  each  radial  line  (at  azimuthal  direc¬ 
tions  i|>  ■  0°,  30°,  60°,  90®,  120°  and  150°)  an  initial  pressure  Pg  so  that  the 
calculated  A-impulse  at  a  large  distance  (e.g.,  at  r  »  30  meters)  matches  with 
the  measured  value.  With  this  set  of  P_  values  extensive  calculation  on  all 
the  physical  quantities  of  interest  are  performed.  The  validity  of  the  BLAST 
code  is  then  judged  from  the  overall  fit  of  the  whole  set  of  calculation  re¬ 
sults  with  the  data. 

The  Pg  values  determined  in  this  way  are  presented  in  Figure  15.  The  ini¬ 
tial  pressurized  balloon  radius  rQ  used  is  again  3  feet.  The  choice  of  rQ  in¬ 
fluences  the  initial  pressure  Pg.  However,  our  calculation  showed  that  the  re¬ 
sults  of  the  calculation,  i.e.,  the  pressure  trace,  A-impulse  and  A-duration, 
etc.,  essentially  do  not  depend  on  this  choice  of  initial  condition. 

Figure  15  also  shows  the  variation  of  Pg  at  the  three  elevation  angles: 
the  shapes  of  the  curves  for  the  three  elevation  angles  are  very  similar,  with 
a  slight  variation  of  the  peak  positions.  This  is  quite  different  from  the 
three  Pg  curves  in  Figure  8  where  the  maximum  pressure  is  used  to  calibrate 
the  initial  condition  of  the  calculation.  The  irregular  variation  of  Pg  in 


Figure  8  reflects  the  effect  of  the  sharp  spikes  which  appear  in  the  pressure 
trace  in  an  unpredictable  fashion.  The  variation  of  Pg  in  the  present  calcula¬ 
tion  (Figure  15),  on  the  other  hand,  reflects  the  change  of  the  flow  field 
near  the  muzzle  brake  when  the  orientation  of  the  muzzle  brake  is  varied. 

In  Figures  16-18  results  of  the  calculated  A-impulse  (solid  curve), 
A-duration  (dashed  curve),  and  maximum  static  overpressure  (dot-dash  curve) 
are  shown  as  functions  of  radial  distance.  The  corresponding  measured  values 
are  represented  there  by  circles,  triangles  and  squares,  respectively.  In 
Figure  19  calculated  pressure-time  histories  are  presented  for  a  few  selected 
cases. 

From  Figures  16-19,  one  sees  that  the  overall  performance  of  the  BLAST 
code  is  good.  The  calculated  A-impulse,  A-duration  and  maximum  static  over¬ 
pressure  follow  the  experimental  values  closely  for  all  azimuthal  orientations 
and  elevation  angles. 

In  comparing  the  calculation  with  the  field  data,  one  should  of  course 
keep  in  mind  that  the  data  contains  a  substantial  amount  of  uncertainty,  as 
Indicated  by  the  large  spreading  of  the  three  sets  of  data  points  in  Fig¬ 
ures  16-18.  The  uncertainty  in  the  data  seems  most  significant  in  the 
A-duration.  This  is  probably  due  to  the  sharp  spikes  which  might  occur  on  the 
descending  slope  of  the  pressure  trace.  We  therefore  feel  that  the  comparison 
for  the  A-duration  is  meaningful  only  in  order  of  magnitude  and  in  this  re¬ 
spect  the  agreement  of  the  calculated  and  experimental  values  is  excellent. 

Figures  16-18  show  that  the  calculated  maximum  static  overpressures  are 
considerably  lower  than  the  corresponding  measured  values,  especially  at 
smaller  distances  (r  *  10  and  20  meters).  On  the  other  hand,  examinations  of 
the  measured  traces  in  Figures  9-11  show  that  a  substantial  part  of  the  mea¬ 
sured  peak  value  at  these  distances  is  contributed  by  the  sharp  spike.  If  one 
discounts  the  spikes  and  compares  the  main  body  of  the  peaks,  excellent  agree¬ 
ment  between  the  calculated  and  measured  maximum  static  overpressure  is  in 
sight.  Comparisons  of  the  calculated  pressure  trace  in  Figures  19a,  b  and  c 
with  the  measured  trace  in  Figures  9c,  10c  and  11c  will  confirm  the  agreement. 

Since  the  calculation  was  performed  so  that  the  calculated  and  measured 
A-impulses  are  matched  at  r  *  30  meters,  comparison  of  this  quantity  should  be 
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concentrated  on  its  values  at  other  distances.  Figures  16-18  show  that  the 
calculated  A-impulse  at  r  ■  10  meters  is  often  smaller  than  the  measured 
values.  We  think  that  this  discrepancy  may  be  due  to  the  nonlinear  effects  of 
ground  reflection  that  are  not  included  in  the  present  treatment.  Since  the 
shock  waves  are  stronger  at  nearer  distances,  the  nonlinear  effects  are  also 
stronger  there.  This  discrepancy  is  not  large,  however,  at  all  and  does  not 
exceed  20Z. 

In  summary,  comparison  of  calculations  and  field  data  presented  in  this 
section  has  shown  that  the  BLAST  code  is  capable  of  describing  the  blast  wave 
generated  by  the  Ml 98  howitzer.  Apart  from  possible  nonlinear  effects  in  the 
near  field,  the  code  satisfactorily  reproduces  the  pressure-time  history, 
maximum  overpressure,  A-impulse  and  A-duration  for  all  azimuthal  orientation 
and  for  the  three  gun  elevations  considered. 

2.3  CATALOG  OF  BLAST  PARAMETERS  AND  WAVE  FORMS 

We  present  here  a  catalog  of  blast  parameters  and  wave  forms  as  obtained 
from  calculations  by  use  of  the  BLAST  code.  This  will  facilitate  quick  refer¬ 
ence  to  the  blast  wave  properties  without  performing  new  computer  runs  with 
the  specific  situation  at  hand.  A  typical  situation  would  be  the  case  where  a 
certain  physical  quantity  such  as  the  A-impulse  or  the  maximum  overpressure  at 
a  certain  location  is  known  and  it  is  desired  to  find  other  information  such 
as  the  pressure-time  history  of  the  blast  wave.  With  the  help  of  the  catalog 
complete  information  can  be  obtained  on  the  physical  quantities  at  that  loca¬ 
tion  and  other  locations  along  the  same  radial  line  from  the  muzzle. 

As  explained  in  Section  2.1,  the  calculation  along  each  radial  line  is  an 
Independent  calculation  on  its  own  right.  The  only  parameters  in  the  calcula¬ 
tion  along  a  radial  line  are  the  height  H  of  the  muzzle  from  the  ground  and 
the  initial  pressure  Pfl.  (The  radius  rQ  of  the  initial  pressurized  balloon  is 
conveniently  chosen  to  be  3  feet  and  the  height  z  of  the  observation  points  is 
fixed  at  5  feet.)  In  Section  2.2  where  the  M198  howitzer  is  considered  speci¬ 
fically,  the  height  of  the  nuzzle  is  determined  by  the  gun  elevation  angle  6. 
In  the  catalog,  it  is  preferable  to  use  the  height  H  instead  of  the  elevation 
6,  since  the  application  of  the  catalog  to  other  guns  is  foreseeable.  Note 
also  that  the  azimuthal  angle  ^  is  not  a  parameter  in  the  calculation.  It  is 


only  a  label  for  a  specific  radial  line  In  the  calculation.  The  catalog 
presents  tables  and  graphs  of  maximum  static  overpressure.  A-impulse  and 
A-duratlon  and  graphs  of  pressure-time  history  for  mizzle  height  H  at  5,  10, 
15  and  20  feet  and  for  Initial  pressures  Ps  at  2,  4,  6,  8,  10  and  12 
atmospheres. 

In  order  to  Illustrate  the  usage  of  the  catalog,  two  examples  are  dis¬ 
cussed. 

Example  1; 

Suppose  in  a  firing  of  a  certain  gun  whose  muzzle  is  10  feet  about  the 
ground,  the  A-impulse  is  measured  to  be  5.0  psi-ms  at  the  100  ft  position  on 
the  line  with  an  azimuthal  angle  of  60°.  With  this  Information,  one  can  locate 
in  the  catalog  the  tables  and  figures  with  the  muzzle  height  value  H  ■  10  ft. 
In  the  table  of  A-impulse  values  (Table  6b)  the  entries  in  the  horizontal  row 
with  r  -  100  ft  can  be  compared  with  the  measured  value  5.0  psi-ms.  In  this 
case  the  value  5.10  psi-ms  is  near  to  the  desired  value  5.0  psi-ms.  Since  this 
value  of  5.10  psi-ms  is  located  in  the  last  column  of  the  table,  it  fixes  the 
initial  pressure  Pg  ■  12  atm.  With  this  Pg  value  one  can  then  find  that  at 
the  same  point  the  maximum  overpressure  is  1.19  psi  (from  Table  6a),  the 
A-duration  is  9.59  ms  (from  Table  6c)  and  the  pressure-time  history  is  shown 
in  Figure  23f.  Information  at  other  distances  along  this  60°  azimuthal  line 
can  also  be  found  under  the  entries  with  Pg  ■  12  atm.  For  example,  at  50  ft 
along  this  line,  the  maximum  overpressure  is  2.50  psi,  A-impulse  is  9.03 
psi-ms,  A-duration  is  7.95  ms,  and  the  pressure-time  history  is  shown  in  Fig¬ 
ure  23f. 

In  the  case  that  the  measure  value  of  A-duration  (e.g.,  5  psi-ms)  does 
not  come  close  to  any  single  calculated  value  but  lies  in  between  the  values 
in  two  adjacent  columns,  one  can  use  interpolation  to  find  a  rough  estimate  of 
the  quantities  discussed  above. 

Example  2: 


The  catalog  can  be  also  used  for  open  charges.  If  the  charge  is  a  spheri¬ 
cal  charge,  the  blast  wave  will  not  depend  on  the  azimuthal  angle  and  the 
values  in  the  catalog  are  applicable  to  all  <|>  values.  In  order  to  make  the 


catalog  more  useful  In  the  case  of  open  charge,  we  give  in  Table  4  the  equiva¬ 
lence  of  pounds  of  TNT  and  the  initial  pressurized  balloon  with  radius  equal 
to  3  ft  and  pressure  equal  to  Pg. 

Suppose  we  have  an  open  charge  of  2  lb  of  TNT  which  is  detonated  15  ft 
above  the  ground.  We  want  to  find  the  blast  wave  properties  at  all  locations 
in  the  field.  Since  1.90  lb  is  quite  close  to  2  lb,  we  use  the  equivalent  P_ 
■  6  atm  for  the  1.90  lb  TNT  for  a  rough  estimate.  (One  can  interpolate  between 
1.90  lb  and  2.65  lb  to  find  the  more  accurate  equivalent  Pg  to  2.0  lb  from 
Table  4.)  Then,  the  maximum  overpressure,  A-lmpulse  and  A-duration  at  various 
radial  distances  from  the  charge  are  given  in  column  3  of  Tables  7a,  b,  and  c, 
respectively.  The  pressure-time  histories  can  be  found  in  Figure  25c. 
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2.47 

3.10 

3.92 

4.82 

120 

1.01 

1.74 

2.31 

2.97 

3.65 

4.45 

* 

c.  A-Duration  (ms) 

•S 

10 

2.26* 

2.41* 

2.59* 

2.77* 

5.10 

4.92 

a 

20 

4.58 

4.77 

5.08 

5.33 

5.63 

5.81 

i 

30 

5.00 

5.16 

5.36 

5.82 

6.06 

6.50 

40 

5.54 

5.67 

5.89 

6.26 

6.59 

6.97 

mV 

‘.1 

50 

6.14 

6.13 

6.38 

6.66 

7.07 

7.40 

60 

6.77 

6.63 

6.78 

7.12 

7.46 

7.89 

N 

70 

7.24 

7.16 

7.37 

7.47 

7.89 

8.27 

V  , 

80 

7.85 

7.70 

7.82 

7.99 

8.20 

8.68 

90 

8.44 

8.22 

8.27 

8.37 

8.67 

8.95 

100 

8.87 

8.76 

8.75 

8.78 

9.02 

9.39 

110 

9.42 

9.28 

9.21 

9.18 

9.37 

9.69 

120 

9.81 

9.65 

9.53 

9.59 

9.73 

10.00 

**In  this  and  the  following  tables,  an  asterisk  denotes  a  case  where  the  pressure 
decreases  to  ambient  before  the  ground  reflection  peak  arrives.  Since  the  A- 

m 

duration 

is  defined  as  the 

time  of  the 

first  passage  of 

the  pressure 

to 

the 

V 

•  • 

ambient. 

the  calculation  of 

the  A-duration  and  A-impulse  counts  only  the  first 

(direct 

incident) 

peak  in  such  a  case. 

Those  values  are  considerably 

smaller 

than  those  where 

both  the  direct  incident  and  the 

ground  reflection  peaks 

are 

;S; 

counted. 

and  are 

not  plotted  in  the  graphs.  The 

maximum 

overpressure 

is 

not 

affected  by  thl8  early  passage  to  ambient. 
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Table  6. 

(H  -  10  ft) 


I 


r 

Pa 

(atm) 

(ft) 

2 

4 

6 

8 

10 

12 

a.  Maximum 

Static  Overpressure 

(psl) 

10 

1.89* 

4.03* 

5.78* 

8.15* 

10.21* 

12.83* 

20 

0.93* 

1.83 

2.44 

3.38 

4.31 

5.13 

30 

0.65 

1.24 

1.78 

2.46 

3.14 

3.75 

40 

0.58 

1.10 

1.51 

1.98 

2.43 

3.06 

50 

0.51 

0.93 

1.28 

1.64 

2.04 

2.50 

60 

0.45 

0.79 

1.07 

1.40 

1.70 

2.04 

70 

0.40 

0.71 

0.94 

1.20 

1.49 

1.77 

80 

0.36 

0.63 

0.82 

1.07 

1.27 

1.54 

90 

0.32 

0.56 

0.73 

0.94 

1.15 

1.34 

100 

0.30 

0.50 

0.67 

0.85 

1.00 

1.19 

110 

0.27 

0.47 

0.61 

0.76 

0.91 

1.07 

120 

0.25 

0.43 

0.56 

0.70 

0.84 

0.98 

b. 

A-Impulse 

(psi-ms) 

10 

1.81* 

4.09* 

6.03* 

8.67* 

10.69* 

13.17* 

20 

1.28* 

4.49 

6.49 

9.72 

12.96 

16.34 

30 

1.83 

3.60 

5.46 

7.69 

10.28 

13.14 

40 

1.67 

3.25 

4.61 

6.37 

8.14 

10.58 

50 

1.50 

2.85 

4.00 

5.48 

6.98 

9.03 

60 

1.43 

2.52 

3.52 

4.81 

6.14 

7.76 

70 

1.35 

2.32 

3.17 

4.22 

5.55 

6.87 

80 

1.27 

2.15 

2.95 

3.94 

4.85 

6.19 

90 

1.21 

2.03 

2.73 

3.59 

4.56 

5.61 

100 

1.12 

1.90 

2.58 

3.34 

4.10 

5.10 

110 

1.05 

1.84 

2.47 

3.09 

3.90 

4.80 

120 

1.02 

1.76 

2.33 

2.95 

3.69 

4.49 

c. 

A-Duration  (ms) 

10 

2.33* 

2.59* 

2.77* 

2.96* 

2.99* 

2.99* 

20 

3.31* 

6.35 

6.36 

6.66 

6.89 

7.02 

30 

6.01 

6.06 

6.33 

6.59 

6.94 

7.35 

40 

6.29 

6.36 

6.54 

6.88 

7.18 

7.55 

50 

6.66 

6.76 

6.84 

7.24 

7.48 

7.95 

60 

7.29 

7.13 

7.26 

7.58 

7.91 

8.33 

70 

7.65 

7.55 

7.60 

7.83 

8.25 

8.62 

80 

8.18 

8.02 

8.13 

8.29 

8.50 

8.97 

90 

8.71 

8.49 

8.54 

8.63 

8.92 

9.34 

100 

9.24 

8.98 

8.96 

8.99 

9.23 

9.59 

110 

9.61 

9.46 

9.40 

9.34 

9.69 

10.01 

120 

10.13 

9.95 

9.84 

9.74 

10.02 

10.30 
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r 

Pg  (atm) 

v'v' 

(ft) 

2 

4  6  8 

10 

12 

a .  Maximum 

Static  Overpi 

41* 

2.94* 

4.31*  5 

85 

1.64* 

2.30*  2 

55 

1.03 

1.43  1 

48 

0.90 

1.24  1 

7.31* 

3.69 

2.39 


2.63* 

3.49* 
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3.  BLAST  LOADING  ON  TWO-DIMENSIONAL  BODIES 


3.1  CALCULATION  OF  GAS  DYNAMICS  ABOUT  IDEALIZED  BODY  SHAPES 

The  coupling  between  the  blast  field  and  the  body  response  hinges  criti¬ 
cally  on  the  forces  Induced  on  the  body  as  the  blast  wave  passes  by.  These 
forces  are  not  constant  around  the  body  because  of  sheltering  effects  due  to 
the  gas  dynamic  flow.  It  is  important  to  any  model  describing  the  motion  of 
the  body  to  be  able  to  describe  the  magnitude  and  timing  of  the  loading  on 
different  parts  for  various  strength  and  orientations  of  the  blast  waves.  In 
this  section  results  are  presented  that  quantify  these  forces  for  use  in  body 
dynamic  calculations. 

The  calculations  were  originally  intended  to  be  compared  with  experimen¬ 
tal  data  taken  at  Lovelace  Research  Institute,  but  because  of  the  inability  of 
Lovelace  to  schedule  tests  during  the  period  of  performance  of  the  contract, 
only  the  calculations  are  presented  at  this  time.  However,  we  have  conducted 
more  calculations  than  originally  called  out  in  order  to  give  a  broader  view 
of  the  possible  loading  distributions. 

In  the  accompanying  figures,  the  transient  loading  of  the  body  due  to  gas 
dynamics  is  presented  for  three  different  body  shapes,  a  square,  a  circle,  and 
a  two-to-one  ellipse  for  two  strengths  of  blast  waves,  3  and  20  psi,  and  for 
two  orientations  of  the  ellipse,  normal  and  45°  rotation.  For  each  geometric 
configuration  and  strength  of  blast  wave,  results  are  presented  for  the  pres¬ 
sure  contours  In  the  air  and  about  the  body,  the  velocity  vector  field  of  the 
flow  generated  by  the  passage  of  the  blast  wave,  and  the  time-dependent  load 
distribution  on  the  body.  The  results  are  self  explanatory  and  will  be  de¬ 
scribed  in  generic  terms. 

In  each  instance,  the  incoming  blast  wave  (indicated  by  the  bunching  of 
pressure  contour  lines  around  the  maximum)  propagates  toward  the  body,  dif¬ 
fracts  around  the  body,  and  creates  a  bow  shockwave  off  the  front  of  the  body 
that  rebounds  toward  the  incident  direction.  The  nature  of  the  orientation  of 
the  body  strongly  affects  the  magnitude  of  the  reflected  pulse,  extreme  exam¬ 
ples  are  between  the  ellipse  at  45°  and  at  normal  direction.  At  normal 


direction  there  is  greater  frontal  area  blocking  the  wave  and  much  more  in¬ 
tense  loading  produced.  The  load  distribution  seen  with  each  geometric  config¬ 
uration  clearly  points  out  the  strong  variation  with  position  and  time  that 
occurs  for  each  body.  In  general,  the  front-facing  part,  which  is  designated 
by  8  *  -90°,  receives  the  largest  impulse  and  is  also  the  first  to  receive  the 
blast  wave.  Approximately  2  to  3  msec  later,  the  time  required  for  the  blast 
to  diffract  around  the  body,  the  back  feels  a  lesser  intensity  loading.  The 
sides  tend  to  see  an  intermediate  value  close  to  the  incident  wave  strength, 
but  one  that  changes  considerably  with  geometry. 

The  shock  multiplying  effect  is  also  seen  in  the  data.  The  lower  level  3 
psi  waves  produce  maximum  frontal  loadings  between  5  and  6  psl  whereas  the  20 
psi  waves  produce  frontal  loadings  as  large  as  70  to  90  psi.  This  amplifica¬ 
tion  is  known  in  strong  shock  wave  theory  and  is  accurately  reproduced  in 
these  calculations. 

The  distributions  presented  here  when  validated  by  experimental  compari¬ 
son,  offer  an  important  link  in  being  able  to  describe  the  body  dynamics.  The 
calculational  tool  used  is  JAYCOR's  EITACC  Code  which  is  capable  of  describing 
arbitrary  geometries,  highly  compressible  flow,  and  uses  a  boundary  condition 
treatment  with  maximum  resolution  and  yet  allowing  waves  to  propagate  out  of 
computational  mesh. 

3.2  PROTOCOL  FOR  VERIFICATION  TESTS  AT  LOVELACE  INHALATION  AND  TOXICOLOGY 

RESEARCH  INSTITUTE  (ITRI) 

This  section  contains  a  "Protocol  for  Exposing  a  Model  of  the  Upper  Torso 
to  Shock  Waves."  The  protocol  was  prepared  by  JAYCOR  after  informal  discus¬ 
sions  with  scientists  at  WRAIR  and  at  Lovelace  ITRI.  At  the  time  of  this  re¬ 
port  preparation,  testing  of  such  a  model  has  been  delayed  due  to  other  and 
higher  priority  tests.  However,  the  basic  objective  and  derived  methodology 
for  the  testing  is  still  valid.  As  a  natural  adjunct  to  the  calculation  of  gas 
dynamics  about  idealized  body  shapes  described  in  Section  3.1,  the  protocol 
should  achieve  the  stated  objectives.  It  is  expected  that  there  may  be  some 
minor  changes  to  the  protocol  before  testing  Is  initiated.  For  instance,  the 
number  of  transducers  to  gather  information  on  the  shock  waves,  the  size  of 
transducers,  method  of  shock  wave  promulgation,  and  others.  These  final  test 
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PROTOCOL  FOR  EXPOSING  A  MODEL 
OF  THE  UPPER  TORSO  TO  SHOCK  WAVES 


1.  OBJECTIVE 

The  objective  of  this  test  procedure  is  to  describe  the  procedures  re¬ 
quired  to  determine  the  effects  of  shock  waves  on  a  model  of  the  human  upper 
torso.  The  requirement  for  testing  is  called  out  in  the  Work  Statement  of  Con¬ 
tract  #DAMD  17-81-C-0104  by  the  Walter  Reed  Army  Institute  of  Research. 

2.  BACKGROUND 

Results  of  the  US  Army  Operational  Test  and  Evaluation  Agency  evaluation 
of  the  M198  howitzer  at  Ft.  Sill  in  July-December  1975  indicated  that  crew  and 
test  personnel  experienced  headaches  and  general  distress  by  being  in  the  im¬ 
mediate  area  of  the  howitzer  when  the  M203  charge  was  used. 

In  1977  the  Surgeon  General  of  the  Army  was  requested  to  examine  the  blast 
effects  of  the  M198  and  other  weapons.  In  the  course  of  that  investigation  the 
US  Army  Aeromedlcal  Research  Laboratory  found  that  the  pressure  levels  exceed¬ 
ed  the  maximum  allowed  by  Mil.  Std.  1474A  for  humans,  and  sheep  exposure  to 
the  shock  tube  with  similar  pressures  at  Lovelace  Research  Institute  showed 
apparent  lung  damage. 

To  better  quantify  the  effects,  detailed  mappings  of  the  weapon  pressure 
fields  were  made  in  May,  1979.  In  October  1979  a  pilot  sheep  study  was  con¬ 
ducted  at  Aberdeen  Proving  Grounds  that  indicated  possible  lung  Injury.  In 
July  and  August  1980  a  detailed  sheep  study  was  carried  out  at  Aberdeen  that 
indicated  acute  lung  and  gastro-intestinal  injury  at  15-17  psi,  but  led  to  no 
definite  conclusion  for  exposures  in  the  2-4  psi  range. 

JAYCOR  has  been  tasked  to  initiate  an  analytical  biomechanical  model  to 
help  guide  and  interpret  exploratory  laboratory  tests  with  animals,  to  esti¬ 
mate  the  damage  risk  criteria  for  humans,  and  to  assist  the  search  for  nonln- 
vasive  procedures.  The  blast  wave  loading  on  a  body  involves  multidimensional 
flow  around  solid  boundaries.  In  order  to  validate  the  results  of  the  computer 
simulation  of  the  blast  interactions  on  a  body,  a  series  of  model  tests  will 
be  conducted  to  determine  the  flow  of  the  blast  wave  about  the  body. 


3.  EQUIPMENT  AND  FACILITIES 


The  torso  model  (Figures  1  and  2)  will  be  fabricated  by  JAYCOR's  San  Diego 
laboratory.  More  discussion  as  to  dimensions  are  discussed  in  Section  III  - 
Pre  Test  Preparation. 

The  facilities,  instrumentation  and  test  personnel  of  the  Lovelace  Inhal¬ 
ation  and  Toxicology  Research  Institute  at  Kirtland  AFB,  Albuquerque,  New  Mex¬ 
ico  will  be  used  to  carry  out  the  tests. 
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SECTION  II 
TEST  PROCEDURES 


4.  GENERAL 

Two  (2)  models  will  be  delivered  to  Lovelace  at  Kirtland  AFB  for  testing. 
Following  mutual  agreement  on  the  selected  transducers,  the  transducer  holes 
and  threading  will  be  accomplished  by  JAYCOR. 

A  series  of  shockwaves  will  be  generated  by  either  bare  charges  or  from  a 
shock  tube  and  allowed  to  impact  against  the  model.  The  model  will  be  rotated 
in  30°  increments  (0°,  30°,  60°,  90°,  120°,  1506,  180°),  see  Figure  3,  with 
surface  pressure  distributions  determined  at  each  orientation. 

The  peak  pressures  of  the  blast  waves  to  be  used  are  25,  10,  5  and  3  psi 
as  measured  with  a  pressure  transducer  not  mounted  on  the  model.  The  pressure 
time  histories  of  the  shots  should  be  similar  to  those  found  in  the  M198  how¬ 
itzer  firings  with  a  fast  rise  time  and  an  A-duratlon  of  approximately  5  ms 
when  in  the  3-10  psi  range.  With  25  psi  it  is  expected  that  tne  A-duration 
will  be  longer.  A  minimum  of  28  shots  (4  pressures  and  7  orientations)  are  re¬ 
quired. 

5.  INSTRUMENTATION 

Seven  transducers  and  their  read  out  instruments  will  be  provided  by  Love¬ 
lace.  The  accuracy  of  the  instruments  should  be  able  to  resolve  the  time-of- 
arrival  differential  among  different  pressure  transducers  on  the  torso  model. 

6.  CALIBRATION 

The  transducers  should  be  calibrated  against  each  other  to  verify  the  in¬ 
formation  on  rise  time  and  sensitivity  provided  by  the  vender.  This  can  be 
achieved  by  exposing  all  the  pressure  transducers,  in  an  array,  to  the  same 
shock  wave. 

7.  EFFECT  OF  SURFACE  MATERIAL 

Mole  skin  is  used  to  wrap  the  model  in  order  to  simulate  human  skin  char¬ 
acteristics  (See  Section  III  Pre-Test  Preparation).  The  pressure  response 
under  the  skin  tissue  is  expected  to  be  different  from  that  on  the  torso 
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surface.  In  order  to  see  such  an  effect,  one  shot  at  a  peak  pressure  of  5  psi 
will  be  tested  before  the  model  is  wrapped  with  the  mole  skin. 

8.  METEOROLOGICAL  DATA 

Before  each  shot  the  ambient  pressure,  temperature,  wind  velocity  and  di¬ 
rection  (if  test  is  conducted  in  open  field)  should  be  recorded. 

9.  SHOT  SCHEDULE 


Round  No 


Peak  Pressure  (psi)  Orientation 


SECTION  III 
PRE-TEST  PREPARATION 


10.  TORSO  MODEL 

The  torso  model  will  be  built  at  the  JAYCOR  San  Diego  Facility.  Although 
only  one  model  is  required  for  the  test,  JAYCOR  will  build  two  identical 
models  so  that  one  will  serve  as  a  spare. 

The  model  will  consist  of  a  torso  cylinder,  two  arm  cylinders  and  two  end 
plates.  The  torso  cylinder  will  be  composed  of  two  semicircular  sections  lo¬ 
cated  at  the  two  sides  adjacent  to  the  arm  cylinders  and  two  flat  sections  at 
the  middle  of  the  front  and  back  side.  (See  Figure  1.)  This  torso  cylinder 
will  be  made  of  aluminum  of  0.5  inch  thick.  The  arm  cylinders  will  be  circular 
in  cross  section.  The  arm  cylinders  and  the  end  plates  will  be  made  of  alumi¬ 
num  or  wood,  whichever  material  is  economical.  The  torso  cylinder  and  the  arm 
cylinders  will  be  wrapped  with  mole  skin  to  simulate  roughly  human  skin  char¬ 
acteristics.  The  end  plates  are  attached  to  the  cylinders  by  screws,  so  that 
they  can  be  easily  taken  apart  for  transducer  Installation  and  services.  A 
hole  (with  a  3  inch  diameter)  In  the  middle  of  the  bottom  end  plate  provides 
the  passage  of  the  connecting  wires  of  the  transducers. 

Seven  transducers  will  be  used  in  measuring  the  pressure  distribution  at 
the  middle  girdle  of  the  torso  cylinder.  Because  of  the  symmetry  of  the  model, 
all  the  transducers  will  be  placed  on  one  side  of  the  torso  cylinder.  (See 
Figure  2.)  The  transducers  and  their  wiring  will  be  provided  by  Lovelace.  To 
our  knowledge,  the  Susquehanna  ST-2  transducer  has  a  pressure  range  1  to  100 
psl  and  is  appropriate  for  the  present  test.  When  this  transducer  is  agreed 
upon  by  Lovelace,  JAYCOR  will  prepare  the  threaded  holes  (hole  diameter:  3/4 
inch;  thread  density:  16  threads/inch)  for  transducer  installation. 
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The  dimensions  of  the  model  are  as  follows: 

h  ■  48  inch 

d^  ■  13.5  inch 

d£  “  10.0  inch 

d3  *  di  -  m  3.5  inch 

d^  “  4.5  inch 

dj  “  0.5  inch 

lj  ■  36  inch 

12  ■  18  inch  only  one  bottom  plate  is  required  for  mounting 

A  table  of  soldier  sizes  from  Military  Standard  1472  is  reproduced  in  Figure 
4.  As  shown,  the  chest  breadth  and  chest  depth  of  the  95  percentile  ground 
troops  are  approximately  equal  to  our  d^  and  d2  values,  respectively.  Our  d^ 
value  gives  an  arm  circumference  of  14  inches  which  matches  closely  to  average 
of  bicep  and  forearm  circumferences  of  the  95th  percentile  ground  troops.  The 
dj  is  chosen  by  rough  estimate. 

The  height  of  the  model,  h,  is  chosen  in  such  a  way  so  that  the  model  is 
free  from  end  effects  and  short  enough  for  easy  handling.  The  dimensions  of 
the  end  plate  are  chosen  by  convenience.  An  error  within  0.5  inch  for  all  the 
dimensions  are  considered  acceptable. 

11.  TEST  SETUP 

Since  the  blast  wave  near  the  end  of  the  shock  tube  is  a  good  approxima¬ 
tion  to  a  plane  wave,  a  test  at  that  location  is  preferable.  In  this  case,  the 
torso  model  should  be  placed  vertically  to  guarantee  normal  incidence  of  the 
blast  wave  on  the  model.  The  model  should  be  supported  by  a  stand,  or  by  other 
means,  so  that  the  center  of  the  model  is  near  the  centroid  of  the  shock  tube 
to  reduce  any  possible  end  effect.  The  orientation  of  the  model  should  be 
changeable  to  facilitate  the  seven  rotations  mentioned  in  Section  II,  Test 
Procedure. 

The  test  with  the  shock  tube  may  be  expensive  and  Lovelace  may  choose  to 

conduct  the  test  with  open  charge.  In  that  case,  the  charge  should  be  suffi¬ 
ciently  far  away  from  the  torso  model  (>  10  d^)  so  that  the  wave  front  at  the 

vicinity  of  the  model  can  be  approximated  by  a  plane  wave.  Assuming  the  model 
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4.  SURVEY  OF  BIOMECHANICAL  MODELING 


4.1  LITERATURE  SEARCH 

In  order  to  properly  evaluate  the  state-of-the-art  in  modeling  of  the 
pulmonary  system,  a  literature  search  was  initiated  to  identify  the  following 
areas: 

•  Markers  for  lung  injury. 

•  Key  scientists  who  have  contributed. 

•  Research  projects  that  might  be  of  assistance  to  Walter  Reed  Army 
Institute  of  Research. 

The  task  was  started  utilizing  the  data  banks  of  NTIS,  BIOSIS  and  NASA.  The 
search  initially  looked  for  lung  and  chest  models  and  associated  experimental 
data.  Of  this  first  search  (Search  I)  we  found  the  following  citations: 


NTIS 

62  citations 

BIOSIS 

137 

NASA 

29 

228 

In  order  to  broaden  the  area  of  the  pulmonary  system,  a  second  search  was 
begun  (Search  II).  Search  II  included  key  words  that  included  the  thorax  and 
injuries  to  this  organ.  This  effort  produced  the  following: 

MEDLINE  158  citations 

NTIS  12 

170 


From  these  literature  searches,  pertinent  articles  were  ordered  and  re¬ 
viewed.  The  relevant  literature  was  then  forwarded  to  WRAIR  for  analysis.  The 
literature  searches  are  Included  as  Appendix  A  and  Appendix  B  to  this  report. 

The  entire  effort  of  four  months  of  literature  search  and  review  of  arti¬ 
cles  was  of  great  assistance  in  identifying  the  state-of-the-art  and  the  major 


contributors.  Many  of  the  contributors  were  invited  and  participated  in  the 
Biomechanical  Workshop  (see  Section  4.2). 

4.2  BIOMECHANICAL  WORKSHOP 

In  preparation  for  the  biomechanical  workshop  a  list  of  potential  con¬ 
tributors  was  assembled  from  the  literature  searches  mentioned  in  4.1.  Using 
the  list  an  analysis  was  made  as  to  what  the  attendees  might  contribute  to  the 
meeting.  From  this  analysis  a  decision  was  made  by  WRAIR  to  explore  the  areas 
of  finite  element  modeling,  fluid  dynamics  applications  to  modeling,  experi¬ 
mental  ta  from  tests  of  frequency  and  intensity  levels  effects  on  the  pulmo¬ 
nary  system,  shock  data  from  animal  testing,  and  spring-dashpot  modeling. 

Site  selection  was  accomplished  prior  to  the  invitation  being  distrib¬ 
uted.  Due  to  the  proximity  of  Albuquerque,  New  Mexico  to  Lovelace  Inhalation 
and  Toxicology  Research  Institute  (ITRI),  it  was  decided  that  Albuquerque 
would  be  well  suited  as  the  workshop  site. 

After  some  exploratory  telephone  calls  were  made  to  candidate  invitees  to 
the  workshop,  an  official  invitation  was  sent  to  those  who  Indicated  a  defi¬ 
nite  interest  in  the  proceedings  and  were  willing  to  present  their  research 
data. 

The  workshop  consisted  of  five  major  presentations  interspersed  with  dis¬ 
cussion  sessions.  The  exchange  of  information  was  lively  and  carried  over  into 
the  evening  hours.  Tape  recordings  and  notes  were  taken  of  the  entire  proceed¬ 
ings  and  a  synopsis  report  prepared.  All  of  these  materials  are  attached. 
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Dear 


JAYCOR  is  pleased  to  extend  you  an  invitation  to  participate  in  our 
BIOMECHANICAL  WORKSHOP  to  be  held  at  the  Classic  Hotel  in  Albuquerque, 
New  Mexico,  December  8-9,  1980.  As  you  have  already  been  contacted  by 
telephone  and  expressed  an  interest  in  participating,  I  would  like  to 
provide  you  some  information  as  to  background,  objective  and 
administration  of  the  Workshop. 

Walter  Reed  Army  Institute  of  Research  (WRAIR)  is  conducting  a 
research  program  in  the  pathophysiology  of  blast  overpressure  in  the 
crew  area  of  military  weapon  systems.  The  possibility  of  nonauditory 
injury  to  soldiers  who  fire  present  or  future  weapon  systems  is  of  major 
concern  to  all.  Exposure  to  overpressures  much  greater  than  those  now 
allowed  for  auditory  safety  may  well  be  harmless.  However,  there  is  no 
positive  verification  of  nonauditory  safety  at  higher  pressure  levels. 
One  of  the  key  areas  that  WRAIR  feels  should  be  addressed  at  the 
Workshop  is:  "Delineation  of  the  mechanisms  of  impact  -  blast  injury 
and  identification  of  the  critical  blast  and  thoracic  parameters  which 
determine  injuries."  The  enclosed  paper  by  '  WRAIR  of  the  Army's 
Technical  Plan,  particularly  Annexes  F,J  and  L  will  provide  a  more 
detailed  insight  into  the  biomechanical  approach.  As  one  of  WRAIR 's 
contractors,  we  are  most  hopeful  that  your  active  participation  in  this 
Workshop  will  assist  us  in  our  efforts  to  model  the  pulmonary  system  and 
to  determine  of  human  response  to  complex  blast  waveforms.  I  have  also 
enclosed  three  papers  by  some  of  the  Workshop  participants,  which  should 
assist  you  to  focus  on  the  objectives  for  our  discussions. 
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We  expect  the  Workshop  to  commence  with  a  "working  breakfast"  on 
Monday  December  8th  and  conclude  the  technical  discussions  about  noon  on 
Tuesday,  December  9th.  For  those  who  care  to  attend,  the  Lovelace 
Biomedical  Laboratory  will  provide  a  tour  and  demonstration  of  their 
large  shock  tube  facilities  at  Kirtland  Air  Force  Base  on  Tuesday 
afternoon. 


We  have  reserved  blocks  of  rooms  at  the  Classic  Hotel  (Telephone 
(505)  881-0000)  for  the  Workshop  participants  for  December  7  thru  9.  As 
some  may  bring  family  members  with  them,  please  notify  the  hotel  of  your 
needs  and  identify  yourself  as  attending  the  JAYCOR  Workshop. 


If  there  are  sufficient  guests  of  the  participants  who  may  desire 
scenic  tours,  we  will  arrange  for  them  a  visit  to  "Old  Town" 
Albuquerque,  the  Sandia  Mountain  Peak  Tramway  overlooking  the  city,  and 
if  possible,  a  trip  to  a  local  Indian  Reservation.  These  arrangements 
are  flexible  and  can  be  varied  to  meet  your  guests'  desires. 
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205  South  Whiting  Street 


Alexandria,  Virginia  22304 


(703)  823-2115 


Please  make  your  own  travel  arrangements  to  Albuquerque.  For  those 
arriving  by  air,  there  will  be  limousine  service  provided  by  the  Classic 
Hotel  from  the  airport. 

JAYCOR  will  reimburse  you  for  your  travel  and  living  expenses  plus 
provide  you  an  honorarium  for  your  participation  in  this  Workshop.  I 
will  discuss  reimbursement  procedures  with  you  during  your  stay.  For 
your  presentation  (15-20  minutes  duration)  we  will  have  available  a 
standard  overhead  vugraph  projector,  a  35  mm  slide  projector  and  a  16  mm 
motion  projector  with  sound. 

If  you  have  any  questions  please  call  me  at  (703)  823-1300, 
Extension  274.  After  December  3rd  you  may  contact  me  through  the 
Classic  Hotel. 


Sincerely  yours. 


Henry  C.  Evans,  Jr. 
Program  Manager 
Fluid  Dynamics  Division 
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XYCOR 


December  7 ,  1 980 


Welcome  ■to  Albuquerque,  the  Classic  Hotel  and  the  Bio¬ 
mechanical  Workshop.  I  hope  that  your  visit  here  -to  "the 
SouthMest  will  be  enjoyable  and  Ihe  Workshop  st imul at ing. 

In  order  to  give  each  of  you  a  chance  to  meet  one  another  in 
a  relaxed  atmosphere,  there  Mill  be  a  "Morking  breakfast"  Monday 
morning  at  8:15  a.m.  in  the  CroMn  Room  on  the  ground  floor.  From 
there  Me  Mill  progress  to  NEWS  II,  our  conference  room,  to 
commence  Mith  the  Workshop.  An  agenda  and  roster  of  our 
participants  have  been  included  in  your  packet. 

Monday  evening  Me  have  planned  a  social  hour  in  the  CroMn 
Room  and  a  dinner  at  the  famous  Maria  Teresa  restaurant  Mhich  is 
located  on  the  edge  of  "Old  ToMn".  If  you  have  a  guest (s)  Mhom 
you  Mould  like  to  invite  to  either  of  these  functions,  please  let 
me  knoM  before  lunch. 

Some  of  you  may  have  a  guest  Mho  may  Mant  to  take  a  tour  of 
Albuquerque  and  its  environs.  Please  let  me  knoM  at  breaktast  so 
that  appropriate  arrangements  can  be  made. 

Once  again,  my  thanks  for  part icipat ing  in  this  project.  If 
I  may  be  of  assistance  during  the  Workshop,  please  let  me  knoM. 


HENRY  C.  EVANS,  Jr. 
Program  Manager 
Fluid  Dynamics  Division 


P.O.  Box  370  -  1401  Camino  Del  Mar 
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December  8,  1980 

Cpt.  (Maj.)  Yancy  Phillips,  M.D.  -  Walter  Reed  Army  Institute  of  Research 

Nonauditory  injury  from  BOP  is  suffered  by  the  air-filled  organs:  sinus, 
inner  ear,  lung,  etc.  The  effects  of  nuclear  explosions  (very  large  BOP)  on 
animals  had  been  described  in  the  1960's  by  Dr.  Richmond's  group  at  Lovelace. 

Based  on  the  nuclear  level  blast  exposure  experience,  the  following  char¬ 
acterization  of  blast  injury  has  been  formed.  The  pathology  of  the  injury  is 
revealed  in:  (1)  hemorrhage  (blood  entering  the  air  passages),  (2)  edema  (in¬ 
creased  water  in  the  lung  and  air  passages),  (3)  emphysema  (enlarging  of  the 
small  air  sacs),  (4)  lacerations,  and  (5)  stripping  of  the  bronchial  epitheli¬ 
um.  The  injury  is  Indicated  pathophysiologically  by  (1)  decreased  pulmonary 
compliance,  (2)  increased  physiological  shunt,  (3)  increased  respiration  rate, 
(4)  decreased  tidal  volume,  (5)  hypoxia,  and  (6)  air  emboli  (air  bubbles  in 
the  blood  stream). 

Mechanisms  that  have  been  put  forward  to  explain  the  coupling  of  the  BOP 
to  the  physiology  include:  (1)  contusion  due  to  chest  wall  acceleration  (most 
of  the  damage  is  beneath  the  chest  wall);  (2)  spalling  effects  (compression 
wave  effects  at  the  boundary  between  dissimilar  materials),  (3)  inertial  ef¬ 
fects  such  as  shearing,  and  (4)  implosion  effects  (local  over-compression  of 
bubbles).  The  external  physical  factors  that  might  be  used  as  indicators  of 
the  dose  strength  are  (1)  peak  pressure,  (2)  duration,  (3)  pressure  impulse, 
(4)  frequency  content,  and  (5)  number  of  exposures. 

In  Lovelace  tests  the  primary  cause  of  death  was  air  emboli  and  a  strong 
dependence  on  the  number  of  exposures  was  shown. 
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In  the  work  of  Clemedson  and  his  colleagues,  increased  lung  weight  was 
used  as  an  indicator  of  damage.  They  found  a  strong  correlation  between  in¬ 
creased  lung  weight  and  chest  wall  velocities  exceeding  15  m/sec  within  150- 
200  psec.  This  corresponds  to  chest  wall  accelerations  in  excess  of  10,000 
g's!  A  second  conclusion  of  this  group  is  that  damage  can  be  produced  by  com¬ 
plex  waves  (arising  from  reflections  in  a  bunker,  say)  and  that  are  only  1/5 
of  the  amplitude  (peak  pressure)  of  classic  blast  waves  required  to  cause 
similar  damage. 

The  BOP  program  will  begin  with  simple  waves,  but  will  be  extended  to  the 
more  important  area  of  complex  waves  caused  by  reflection  in  bunkers  or  vehi¬ 
cles. 

The  present  design  criteria  is  contained  in  HILSTD  1474.  The  parameters  of 
this  standard  are  peak  pressure  and  B-duration  (the  time  from  the  start  of  the 
wave  until  the  amplitude  is  20  dB  of  the  peak  for  the  last  time). 
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Below  the  W-line  no  hearing  protection  is  required.  Various  amounts  of 
hearing  protection  is  required  up  to  the  Z-line.  No  exposure  above  the  Z-line 
is  permitted  because  of  possible  nonauditory  damage. 

[Von  Gierke,  who  was  on  the  committee  that  drafted  the  standard,  pointed 
out  that  the  Z-line  was  only  a  conservative  guess.] 

The  variance  from  MIL-STD-1474  by  new  weapons  has  a  serious  impact  on  the 
Army's  function  and  thus  the  urgency  of  the  BOP  program.  The  M198  will  be  the 
principal  field  piece  of  the  1980-90's.  Its  muzzle  brake  redirects  the  blast 
toward  the  crew  area  and  triples  the  BOP  there.  There  is  currently  a  12-month 
moratorium  on  crew  training  using  the  high  range  charge.  A  self-propelled  how¬ 
itzer  with  a  more  efficient  muzzle  brake  (20%  vs.  10%)  has  shown  Z-line  cross¬ 
ings  with  zone-7  charges.  [Von  Gierke:  data  from  other  weapons  should  be 
collected.]  [Cummings:  other  weapons  that  may  exceed  the  Z-line  are  4.2”  mor¬ 
tar,  81mm  mortar.]  WRA1R  will  undertake  field  demographic  data  -  a  cross  sec¬ 
tion  prevalence  study  of  pulmonary  function  in  active  duty  artillerymen. 

Dr.  David  Viano  -  General  Motors,  Research  Laboratory 

General  Motors'  concern  is  with  the  impact  loading  on  the  human  chest  - 
steering  wheel  impact  during  car  accidents. 

Principal  models  have  lumped  mass  and  spring  concept. 
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Fig.  I.  Mechanical  analog  for  the  anteroposterior  impact  re¬ 
sponse  of  the  human  thorax.  Lobdell  [1]. 


The  main  test  of  the  model  is  the  prediction  of  the  sternum-backbone  displace¬ 
ment.  The  impact  leads  to  chest  wall  compression,  whole  body  motion,  and  ener¬ 
gy  dissipation.  No  correlation  of  AIS  with  impact  momentum  was  found,  but 
there  was  correlation  with  kinetic  energy.  It  should  be  noted  that  AIS  is  ex¬ 
tremely  crude  and  subjective  and  is  directed  toward  severe  crushing  injuries. 


Fig.  13.  Dependence  of  cadaver  injury  data  on  the  impactor 
kinetic  energy. 


Fig.  14.  Linear  dependence  of  resultant  injury  on  peak  nor¬ 
malized  deflection  from  cadaver  tests. 

When  rib  damage  is  removed  from  consideration  the  results  indicated  that 
visceral  damage  does  not  occur  until  the  chest  is  40%  compressed.  Compression 
also  serves  to  regroup  the  velocity  correlation. 
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There  appear  to  be  two  injury  modes:  crushing  impact  (relatively  slow  but 
large  displacements)  and  impulse  impact  (small  but  quick  displacements). 


Experiments  were  conducted 

ity. 


t 

Two  classes  of  lung  injury  were  identified:  (1)  alveolar  and  (2)  bronchial. 
Alveolar  injury  is  associated  with  impulse  impact;  bronchial  injury  with 
crushing  impact.  Bronchial  injury  related  to  damage  at  the  root  of  the  lung. 
Alveolar  injury  is  morphologically  similar  to  blast  injury. 

Results  of  experiments  with  constant  velocity  for  a  fixed  distance. 

Latest  modeling  effort  is  directed  toward  a  more  fundamental  understanding 
using  a  finite  element  dynamic  model. 
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to  control  both  impact  displacement  and  veloc- 
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Present  model  is  planar  (two-dimensional)  and  treats  the  rib  cage  as  a  contin¬ 
uous  loop.  Some  difficulties  with  present  model  for  large  geometric  displace¬ 
ments,  so  GM  is  sponsoring  further  work  with  University  of  California,  Berke- 


Dr.  H.  E.  Von  Gierke  -  Air  Force  Aerospace  Medical  Research  Laboratory 


The  Air  Force  Is  concerned  with  Impact  to  pilots  during  ejection  from  jet 
planes.  Tests  have  been  conducted  using  controlled  frequency  vibrations  to  hu¬ 
mans  in  a  sitting  position.  A  modulated  air  stream  from  the  mouth  is  observed. 
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The  same  experience  has  been  observed  with  air  transmitted  forces:  infrasound 
corresponds  to  vibration,  blast  waves  to  impulse. 
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Possible  important  difference  between  Impact  and  blast  loads:  impact  will 
couple  with  low  frequency  modes,  blast  with  high  frequency  modes. 
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Compressibility  of  the  lung  had  to  be  included  in  order  to  explain  the  higher 
frequency  resonance.  The  resonance  at  4-6  Hz  seen  earlier  in  vibration  tests 
does  not  appear  here  because  the  distribution  of  the  load  is  different. 

General  conclusions  on  the  nature  of  response  and  damage : 
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Looks  similar  to  Lovelace  lethality  curves. 

(2)  General  (nonlocal)  damage  will  result  from  the  low  frequency  part  of  the 
wave  whereas  local  damage  will  result  from  high  frequency  part. 

Dr.  Ints  Kalep  -  Air  Force  Aerospace  Medical  Research  Laboratory 

AF  developed  a  detailed  lumped  mass  thorax  model  to  deal  with  problem  of 
spine  compression  during  pilot  ejection.  Properties  of  the  model: 

(1)  lumped  parameter  torso  model  for  Z-vibration  and  impact, 

(2)  transverse  coupling  of  chest  wall, 

(3)  nonlinear  compressibility  and  air  passage  resistance, 

(4)  effects  of  blunt  impact  on  chest  wall. 

Conclusions  reached  were  that:  for  short-duration  impacts  air  dynamics  are  not 
Important  and  that  rib  fractures  correlate  best  with  peak  chest  deflection. 
For  large  chest  wall  deflections  the  resistance  of  the  lung  pressure  must  be 
Included  to  get  the  observed  effects. 
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The  Air  Force  feels  that  certain  model  improvements  are  necessary 
partmentalized  thoracic  volume,  better  tracheal  flow  characteristic! 
effective  chest  wall  elasticity  for  various  pressure  distributions. 


December  9,  1980 

Dr.  Paul  Chen  -  TRW  Defense  and  Space  Systems  Group 

Dynamics  of  the  human  thorax.  The  physiological  components  to  be  modeled 
are  superficial  tissues,  ligate  muscles,  bone  structure,  and  internal  organs. 
The  objectives  of  TRW's  interest  are  crash  impact  injury,  auto  restraint  sys¬ 
tems,  and  antropological  dummy  development  and  manufacturing.  The  problems 
associated  with  human  thorax  modeling  are  injury  criteria,  required  complex¬ 
ity,  and  costs.  The  approaches  considered  by  TRW  are:  (1)  pathophysiological  - 
using  clinical  research  to  determine  the  injury  mechanism  and  a  damage  crite¬ 
rion;  and  (2)  biomechanical  -  using  a  sequence  of  animal,  cadaver,  human  vol¬ 
unteers,  and  anthropomorphic  dummy  testing  coupled  with  analysis  and  inter¬ 
pretation  using  a  mathematical  model.  The  available  modeling  techniques  are: 
(1)  simplified  lumped  parameter  mot 2ls  -  easy  to  construct,  inexpensive,  model 
parameters  are  critical,  and  there  is  usually  limited  information  available; 
and  (2)  detailed  finite  element  model  -  more  complex,  cost  more,  some  model 
parameters  may  still  be  unknown,  generates  detailed  predictive  information. 

Finite  Element  Model.  The  input  quantities  required  are:  body  geometry, 
material  properties,  stiffness,  joint  boundary,  load  distribution,  and  the 
energy  dissipation  coefficients.  The  model  implications  include:  a  skeletal 
module,  viscera  inclusion,  physiological  effects  of  intrathoracic  pressure  and 
muscle  tension.  A  synthesis  technique  is  used  where  a  substructure  of  models 
are  solved  individually  (nonlinear  effects  ignored)  and  then  the  total  solu¬ 
tion  constructed  from  a  combination  of  modes.  The  existing  finite  element  com¬ 
puter  codes  are  ADINA,  SAP,  and  NASTRAN,  but  more  specific  codes  need  to  be 
developed.  The  current  finite  element  chest  model  used  by  TRW  was  originated 
at  UCLA  by  Chen  and  Roberts  with  collaboration  with  Raddi  and  Kazemieslamia. 
THORAX-I :  static,  elastic  model,  high  resolution  (80  nodes?).  THORAX-II :  sim¬ 
plified  model  (~  20  nodes). 

Suggestions:  (1)  Injury  criteria  and  injury  mechanism  study  should  be 
done  in  parallel;  (2)  use  detailed  finite  element  model  that  can  then  be  sim¬ 
plified  to  determine  the  parameters  of  simpler  lumped  models;  (3)  animal 
models  tested  first,  then  human  models;  (4)  a  two-level  model  should  be  used 
with  a  bone  thorax  part  and  a  soft  tissue  part;  (5)  first  use  a  linear  model, 
then  include  nonlinear  effects. 

201 


4.3  LUMPED  PARAMETER  MODELS 


As  part  of  our  investigation  of  the  state  of  biomechanical  modeling,  we 
implemented  existing  lumped  parameter  models  being  used  to  describe  thorax  and 
abdominal  motion  into  a  computer  code.  In  particular,  we  wrote  a  computer  pro¬ 
gram  to  solve  and  graphically  display  the  results  of  the  Lovelace  Model  that 
has  been  used  in  nuclear-level  blast  interaction.  This  computer  program  was 
delivered  to  Walter  Reed  and  put  on  their  in-house  machine  to  allow  more  rapid 
turnaround  on  answering  questions  relating  to  blast  interactions.  We  also  have 
maintained  the  computer  program  on  the  JAYCOR  computer  to  be  able  to  answer 
questions  relative  to  possible  body  motion  under  loading. 

One  of  the  applications  that  has  been  made  of  the  code  is  to  compare  the 
effects  of  a  blast  signature  as  measured  in  the  field  from  a  howitzer  with  the 
loading  signature  developed  by  the  water  jet  impactor  we  are  developing  for 
Walter  Reed.  As  an  example,  there  was  interest  in  knowing  the  effects  of  an 
after  pressure  following  the  initial  pulse.  This  after  pressure  could  be  due 
to  the  mechanical  aspects  of  the  impactor  or  to  the  winds  that  follow  a  blast 
wave.  The  accompanying  figures  show  a  comparison  between  the  response  of  the 
lumped  parameter  model  scaled  to  the  mass  of  a  man  for  the  cases  when  the  in¬ 
cident  blast  wave  does  or  does  not  have  an  after  pressure  part.  In  all  cases, 
the  after  pressure  was  assumed  to  be  10%  of  the  maximum  peak  pressure*  The  re¬ 
sults  are  shown  for  3,  10,  25,  and  50  psi  maximum  peak  pressure  waves.  The 
results  indicate  that  the  responses  are  nearly  identical  and  that  for  present¬ 
ly  accepted  injury  indicating  quantity,  chest  wall  velocity,  there  is  only  a  2 
to  4%  variation  caused  by  the  presence  of  an  after  pressure.  This  variation  is 
so  small  that  it  is  unlikely  that  the  after  pressure  has  any  influence  on 
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5.  DISCRETE  REPRESENTATION  OF  BODY 


5.1  INTRODUCTION 

The  Finite  Element  Method  (FEM)  has  been  used  extensively  in  different 
physical  and  engineering  fields.  For  problems  involving  irregular  geometry, 
boundary  conditions,  prescribed  loadings,  as  well  as  complex  material  proper¬ 
ties  FEM  analysis  furnishes  approximate  solutions  to  the  physical  field  prob¬ 
lems  with  satisfactory  accuracy. 

We  have  begun  to  study  the  lung  injury  mechanism  during  airblast  over¬ 
pressure  by  the  FEM  technique.  By  assuming  linear  material  with  small  de¬ 
formation,  a  two-dimensional  plane  strain  discretization  model  accounting  for 
the  composite  materials  and  complex  structure  of  human  thorax  cross  section  is 
used.  At  this  stage,  static  analysis  has  been  made  to  study  the  deformation 
and  stress  distributions  of  the  structure  model  when  a  static  pressure  (50  psl 
or  5  psi)  is  applied  on  the  frontal  chest  wall,  side  wall,  or  back  wall  of  the 
thorax.  Transient  analysis  of  the  model  response  during  and  following  5  msec 
of  step  pressure  loading  (5  psi)  on  the  front  chest  wall  is  also  being 
studied. 

The  purpose  of  the  current  task  is  the  construction  of  a  plane  strain 
two-dimensional  FEM  model  and  the  investigation  of  structure  deformations 
(static  and  transient)  resulting  from  prescribed  external  pressure  loadings. 
This  is  only  the  initial  phase  of  the  whole  study.  The  objective  will  be  the 
construction  of  a  continuum  model  which  not  only  confirms  the  measurable  ex¬ 
perimental  data  from  various  types  of  field  study  but  also  makes  clear  def¬ 
inite  risk  criteria  judgement.  To  accomplish  this  the  propagation  of  blast 
waves  in  different  parts  of  the  body,  the  scattering  and  diffraction  of  the 
elastic  waves  through  different  media,  and  how  tissue  is  damaged  when  the 
resultant  stress  reaches  the  ultimate  strength  are  to  be  studied  with  the  aid 
of  the  FEM  model. 
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DESCRIPTION 


The  assumptions  used  in  this  analysis  are  summarized  as  follows: 

(1)  2D  Model  of  the  Thorax  -  Based  on  the  cross  section  of  a  human  trunk 
at  the  level  of  the  aortic  semilunar  valve  a  2D  plane  strain  model  of  the 
thorax  is  constructed.  Physiologically,  the  pleural  space  will  allow  relative¬ 
ly  free  sliding  motion  between  lungs  and  chest  wall,  diaphragm,  as  well  as 
other  organs.  For  simplicity,  it  is  assumed  that  there  is  no  relative  motion 
between  surfaces  of  different  organs  during  any  given  loading  and  deformation 
process.  In  other  words,  a  displacement  compatible  model  is  used.  The  FEM  mesh 
discretization  is  shown  in  Figure  37(a)  and  37(b).  The  rib  cage  is  modeled  as 
a  closed  ring  along  with  costal  cartilage  and  sternum  in  model  A.  In  model  B, 
the  rib  is  modeled  as  segmental  with  skeletal  muscle  filling  the  space  between 
to  account  for  the  average  properties  and  roles  of  a  rib  cage  structure. 

(2)  Linear  Analysis  -  It  is  a  well  known  fact  that  most  of  the  biologi¬ 
cal  soft  tissues  are  nonlinear  viscoelastic.  For  nonlinear  material  under 
finite  deformation  the  analysis  is  extremely  involved  and  complicated.  While  a 
nonlinear  three-dimensional  analysis  is  possible,  a  linear  approach  is  used  at 
the  present  stage  for  the  sake  of  simplicity. 

(3)  Material  Properties  -  For  this  linear  analysis,  Hookean  type  materi¬ 
al  constants  are  used.  For  the  lung  parenchyma  and  the  aortic  vessel  wall  they 
are  taken  from  Radford  (1957)  and  Bergel  (1972),  respectively.  For  the  rest  of 
the  materials  the  moduli  constants  are  approximated  from  stress-strain  curves 
collected  in  Yamada  (1970).  G,  K,  E  and  v  denote  shear  moduli,  bulk  moduli. 
Young's  moduli  and  Poisson  ratios,  respectively.  The  magnitudes  of  densities 
are  chosen  arbitrarily  in  this  analysis  with  2.67  x  10  Newton-sec  /m  for  all 
the  tissues  except  the  lung  where  20Z  of  the  above  number  is  used  since  the 
purpose  is  to  study  the  role  of  densities  difference  in  inertia  terms.  Table  9 
summarizes  the  material  property  constants  used. 

(4)  FEAP  Code  -  The  FEAP  code  (Finite  Element  Analysis  Program)  devel¬ 
oped  by  R.  L.  Taylor  of  the  University  of  California  at  Berkeley  is  used  to 
perform  the  present  analysis.  FEAP  is  a  versatile  general  purpose  FEM  program 
with  emphasis  on  its  capability  on  contact-impact  problems. 
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(a)  Model  A  -  Rib  Cage  Modeled  as  a  Closed  Ring  with  Costal  Cartilage 
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(b)  Model  B  -  Rib  Cage  Modeled  as  Segmental  with  Muscle  In  Between 

Figure  37.  TWo  Dimensional  Plane  Strain  FEM  Mesh  Discretization  Shown 
with  Different  Material  Zones* 


Figure  37(c).  Cases  S-3,  S-4,  and  S-5  are  Three  Different  Loading  Cases  on 
Model  A  as  Indicated.  Case  S-6  is  the  case  with  5  psl  of 
front  chest  loading  on  Model  B. 


Table  9.  Material  Property  Constants 


Shear 

Bulk 

Young ' s 

Poisson 

Density 

Modulus, 

Modulus , 

Modulus, 

Ratio, 

Estimated 

G 

K 

E 

V 

P 

Tissues  (Newton/m2) 

(Newton/m2) 

(Newton/®2) 

(Newton-sec2 /m4 ) 

Skeletal  muscle  1.1  *  105  1.1  x  io6  3.3  x  io5  0.449  2.67  x  io3 


Rib  (bone) 


5.0  x  10' 


ao 


1.1  X  I0iv  1.3  X  10 


10 


0.301 


Heart 


7.2  x  10" 


7.0  x  10- 


2.1  x  10- 


0.450 


Blood 

Lung 

Aortic  wall 
Mediastinum 


1.4  x  10' 


1.3  x  10J  3.1  x  10- 


2.8  x  io5 
■»4 


2.8  x  10°  8.0  x  10- 


2.6  x  10-  3.1  x  10J  7.6  x  10" 


0.105 

0.453 

0.459 


2.67  x  ioJ 


Costal  cartilage  1.6  x  108  7.5  x  108  4.5  x  108  0.400  2.67  x  103 


2.67  x  io- 


8.3  x  104  4.2  x  107  2.5  x  105  0.499  2.67  x  103 


5.34  x  io2 


2.67  x  io- 


2.67  x  io3 


In  static  analysis  the  standard  Gauss  elimination  technique  is  used  to 
solve  the  force-deformation  relationship  with  a  triangular  decomposition  of 
stiffness  matrix  K.  In  transient  analysis  direction  integration  scheme  with 
Implicit  solution  is  followed  by  using  a  one-step  Newmark  method  to  discretize 
in  time  and  a  Newton  method  to  solve  the  problem. 

With  FEAP  extensions  of  the  present  study  into  cases  with  nonlinear  elas¬ 
tic  or  linear  viscoelastic  materials  are  possible. 

RESULTS 

For  static  analysis,  three  loading  cases  were  performed  on  model  A. 
Pressure  loading  of  50  psi  is  applied  on  the  frontal  chest  wall,  side  wall,  or 
posterior  wall  in  case  S-3,  S-4,  or  S-5,  respectively. 

On  model  B  (see  Figure  37b),  instead  of  a  closed  ring  type  rib  cage  a 
segmental  type  rib  is  used  to  account  for  the  average  roles  and  properties  of 
a  three-dimensional  rib  cage  structure.  Five  psi  of  pressure  loading  is  ap¬ 
plied  on  the  frontal  chest  wall  (case  S-6)  and  side  wall  (case  S-7).  Load  con¬ 
figurations  are  shown  in  Figure  37c.  Undeformed  as  well  as  deformed  cross 
sectional  geometric  configurations  for  cases  S-3,  -4,  -5,  and  -7  are  shown 
in  Figures  38,  39,  40,  41,  and  42,  respectively. 

For  case  S-3  the  magnitudes  of  the  average  stress  components  at  different 
element  in  each  of  the  organs  are  shaded  in  different  darkness  (Figure  43). 
Follow  the  order  of  lung,  skeletal  muscle,  ribs,  costal  cartilage,  heart, 
aortic  vessel  wall,  mediastinum,  and  blood,  stress  components  ox,  oy,  and  o^ 
are  labeled  on  different  graphs. 


DISCUSSION 

Results  from  static  analysis  give  the  developed  stress  due  to  deformation 
as  the  pressure  loading  is  applied  on  the  external  surface.  Regional  differ¬ 
ences  in  different  organs  in  various  loading  cases  can  be  seen. 


The  magnitude  of  the  stress  in  the  lung  predicted  based  on  model  A  is 
only  of  the  order  0.1  psi  when  the  externally  applied  pressure  is  50  psi, 
since  the  high  stiffness  closed  ring  type  ribs  in  model  A  essentially  coneti- 


(a)  Model  A  with  Loading  Indicated 


(b)  Undeformed  Configuration 


(c)  Deformed  Configuration 
Figure  38*  Case  S-3  In  Various  Configurations. 


50  psl 


(c)  Deformed  Configuration 


Figure  39.  Case  S-4  In  Various  Configurations 


(c)  Deformed  Configuration 


Figure  AO.  Case  S-5  in  Various  Configurations 


(a)  Model  B  with  Loading  Indicated 


(b)  Undeformed  Configuration 


(c)  Deformed  Configuration 


Figure  41.  Case  S-6  In  Various  Configurations 


Figure  42.  Case  S-7  in  Various  Configurations 


organs  Inside  it.  To  account  for  the  roles  of  a  three-dimensional  rib  cage 
during  pressure  loading  we,  therefore,  use  a  segmental  type  rib  model  B  in 
this  two-dimensional  model.  The  idea  is  so  that  average  properties  and  roles 
of  the  bony  and  muscle  parts  can  be  estimated.  With  model  B  the  stress  in  the 
lung  is  of  the  order  of  0.5  psi  when  the  external  pressure  loading  is  5  psi. 
Based  on  this  comparison  it  is  reasonable  to  assert  that  model  B  is  a  better 
candidate*  However,  further  development  and  modifications  of  the  two- 
dimensional  model  is  still  under  way.  We  have  observed  the  stiffness  differ¬ 
ence  of  the  whole  thorax  structure  is  tremendous  as  we  change  from  model  A  to 
model  B.  The  relative  orientation  of  ribs  and  the  lung  tissue  behind  them,  the 
contributing  factor  of  the  rib  cage  to  the  lung  injury  mechanism,  protective 
or  negative?  These  are  all  interesting  topics  to  be  studied. 

We  have  shown  one  way  to  estimate  the  stress  distributions  in  the  lung 
due  to  gross  deformation  resulting  from  pressure  loading  on  the  external  body 
surface.  The  order  of  magnitude,  however,  is  considerably  small.  Experimental 
results  reported  by  Bowen  et  al.  (1968)  indicate  the  intrathoracic  overpres¬ 
sure  of  a  10  kg  dog  could  go  as  high  as  hundreds  of  psi  during  exposure  to  air 
blast  shock  wave.  This  type  of  overpressure  or  high  stress  is  probably  caused 
primarily  by  the  propagation  of  incident  waves  and  their  scatterings  and  dif¬ 
fractions.  Wave  study  is  currently  under  way. 

We  have  considered  the  lung  as  an  organ  made  of  very  compliant  spongy 
material.  In  actuality  the  lung  is  a  collection  of  millions  of  tiny  alveoli 
(air-containing  sacs)  at  a  more  refined  scale.  As  an  air-containing  structure 
during  rapid  compression  the  pressure  volume  relationship  of  the  air  contained 
in  the  lung  should  change  drastically.  The  incorporation  of  a  pressure-volume 
relationship  for  the  air  inside  the  lung  into  the  FEM  model  is  currently  being 
done. 

The  current  model  has  assumed  no  relative  motion  among  interfaces  of  the 
organs.  This  assumption  will  influence  the  loading-response  prediction  and  the 
thorax  model  in  being  extended  to  include  slip. 

We  have  also  performed  transient  analyses.  Sample  problem  of  uniform 
stretching  a  metal  sheet  with  a  center  hole  has  been  carried  out  as  a  test  of 
dynamic  analysis  capability  of  the  FEAP  code.  Transient  analysis  of  our  thorax 


model  B  response  to  5  msec  of  pressure  step  loading  (5  psl)  has  also  been  car¬ 
ried  out.  The  task  of  detailed  transient  response  analysis  is  In  progress. 
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5.2  FINITE  DIFFERENCE  REPRESENTATION  OF  PRESSURE  WAVES  IN  THE  PARENCHYMA 

In  order  to  evaluate  the  effects  of  pressure  waves  in  the  spongy  lung 
parenchyma,  the  JAYCOR  SPUNG  Code  was  used  to  calculate  the  refraction  and 
reflection  of  pressure  waves  in  an  idealized  lung  cross  section.  The  code 
solves  the  material  properties  describing  compressive  wave  propagation  in 
materials  of  varying  properties,  and  was  set  up  to  describe  lung  shapes  re¬ 
gions  of  high  compressibility,  surrounded  by  material  of  low  compressibility 
representing  the  solid  and  liquid-filled  organs  of  the  thorax  cross  section. 
The  model  was  then  excited  by  pressure  waves  of  varying  strengths  passing 
around  the  perimeter  and  observing  the  body  response. 

The  accompanying  figures  show  examples  of  the  complex  pressure  distribu¬ 
tions  that  develop  within  the  lungs  and  the  pressure  time  histories  at  partic¬ 
ular  points  within  the  lung.  The  results  show  a  complex  overpressure  field 
which  locally  can  become  many  times  greater  than  the  incident  waves,  and  can 
lead  to  local  damage,  especially  upon  reflection  from  dissimilar  materials. 
Such  behavior  is  possibly  the  source  of  localized  damage  that  is  thought  to 


Figure  44.  Calculation  of  pressure  wave  propagation  through  a  two-dimensional  cross-section 
of  a  heterogeneous  body  with  stiff,  high  density  sections,  representing  muscle  and 
bone;  and  spongy,  low  density  sections  representing  the  lungs.  The  sequence  of 
pictures  shows  a  wave  incident  from  the  top  left  that  leads  to  an  overpressure 
focussing  in  the  right  lung  that  has  three  times  the  pressure  of  the  incident  wave. 


MAXIMUM  LEFT  LUNG 


Time  History  of  Maximum  Overpressure  In  the  Lung  Modeled  by  the 
SPUNG  Code.  The  Initial  blast  wave  has  a  maximum  pressure  of 
3  psl  yet  local  overpressures  of  almost  10  psi  can  be  formed  by 
wave  focussing.  This  enhancement  of  pressure  and  the  period  of 
oscillation  arise  out  of  the  assumed  material  properties,  which 
In  this  case  are  somewhat  arbitrary  but  could  be  determined  more 
precisely  by  controlled  laboratory  measurements.  There  is,  how¬ 
ever,  qualitative  agreement  with  the  measured  traces  reported 
by  Lovelace. 


6.  SUMMARY  AND  CONCLUSIONS 


In  summary,  the  tasks  performed  under  this  contract  and  its  modifications 
have  scoped  and  developed  the  technologies  required  to  assist  the  blast  over¬ 
pressure  program  in  quantifying  and  validating  the  damage  risk  criteria  for 
future  weapon  deployment*  As  was  mentioned  earlier,  this  effort  has  been 
directed  toward  guiding  the  use  of  technology  as  well  as  implementing  it  into 
a  specific  form.  The  availability  of  complete  and  consistent  models  connecting 
the  blast  source  through  the  far  field  propagation  to  the  loading  on  the  body 
to  the  structural  response,  and  finally  to  the  local  stress  distributions 
within  the  lung  is  an  important  part  of  the  blast  overpressure  program 
mission. 
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INPFM  II :  A  Computer  Implementation  of  Recent  Models  for 

Estimating  the  Oose  Equivalent  to  Ot gano  of  Man  from  an 
Inhaled  or  Itigested  Radionuclide 


Abstract:  Single  detector  arrangements  are  used  to  measure 
be.-v  t  radioactivity  curves  In  healthy  subjects  and  in  patients 
with  various  heart  failures.  A  method  Is  developed  from  a 
modified  gamma  function  to  det*>rmtre  the  cat  d  I  opu  I  monai  y 
parameters  from  the  radlocatdlograms :  systemic  flow,  pulmonary 
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SWcOISf  NULL  h  f  K  PLANT  AT  EARS  EE  ACK.  YhE  RELEASE  FRACTIONS  FOR 
Tnc$f  ACCIDENTS  ARE  TAKEN  FROM  WASh-1430.  BASED  ON  PARAMtTRIC 
STUD  I  L  S*  DOSES  ARE  CALCULATED  FOR  VtR  Y  UNF  AVI  URABL  h*  PUT  NOT 
lNC.RcDl.iLt  WtAlHtR  CONDITIONS.  Tnt  PROBABILITY  OF  SUCn  CONDITIONS 
IN  COMBINATION  WITh  WIND  DIRECTION  TOWARDS  QANlSh  TERRITORY  IS 
t ST  IM  A  T tU «  DOSES  TO  RONE  MARROW*  LuNGS*  G I  — T R A CT  AND  ThYROID  ARc 
C  aLC  JL  fTcU  USING  DCSt  MCl'tLS  OtVtLOPtD  AT  RIS9.  TutSE  DOSES  ARE 
FOUND  ?0  CONSISTENT  WITn  DOSES  CALCULATED  WITh  ThE  MCDFlS  USED 
IN  WAln-l^CO. 

*aCCIl  ’:NTS,  DCCUPATII  NAL  ;  »A1R  POLLUTION*  R AD  10 AC TI  Vfc /ANALY  S  1  S 
BCNc  maRROw/KALI AT  ION  tFPfcCTS  :  COm°UTCKS  ;  DENMARK 

ENVIRONMENTAL  EXPOSURE  ?  GaSTRUINTFST INAL  SY ST tM/K AD  I  AT  ION  EFFtCTS 
HUMAN  ;  LUNG/RA.D1AT1CN  EFFECTS  5  M?Th  tMAT  iCS  J  MCDtLS  *  TntCRtT  1CA  L 
*NUCLIAK  REACTORS  ;  RADIATION  DOS AGt  ;  RISK  ;  SWEDEN 

Thyroid  gland/radi  at  ion  effects  ;  «eathpr 

ENG 

RISC  KcP  1«*77  OCT;  OSM  tl-i9 


0  FPL  INS 


PRINT 


n 

1  AU 

T  1 

'n 

*  % 

At1. 

Mh 

Mh 

^  • 

S' 

Mh 

£ 

mh 

Mh 

la 

SO 

kN! 

2  AU 

i 

T  I 

■ 

A  B 

CASc« 
■SlaST 
and  c 
f  ive 

txPuC 

Wc<t 

Tnc  F 

PaTIE 

V  t\T  1 

D  I  tD  i 

A  N  A  E  >, 

HYPO* 

FuNCT 

A«-TEK 

adult 

CARBU 

LU"''G/ 

J  \t  ►:  jj  v» 

=>  c:  s  1  T 

t\'G 

IvjuR 

Ci’P°t 
y  cast 
up  UN 

CL“i-wiu 

P  fl  t  i  c 

BULL? 

UNMUL 

tXPtP. 

I  b  K  A 

T  HO  US 

BLAST 

HUMAN 

\L‘K  1  n 

•.’LUND 

eNG 

PR  J 


*  vs  ; 

INJUR 
l  UR  S  t  • 
PaT  I  E\t 
S  lUX'S 
ApF>ARt 

x«>losi 

NTS  Wc 
LAT lUN 
CNF  C 
7  hE  S  I A 
*cMIA 

ICi  IE 
TmE  B 
;  AINU 

\  diox 

* IN JUK 
07  HOP  A 
ivt  PR 


PCRTrR  M F 

1Tb  7t'  lt-t  LUNGS:  CLINll.AL  PRt  S  tN  T  Al  1  ON  »  MANAGtMtNl 


TS  WITh  b L 
ARE  RtPOR  7 
\'7  UN  »hc 
PNS.  AK7FK 
K  E  Adi  VtL 
t  WhICh  wA 
WING  10  F  1 
t  AND  7 1-iE 
PERSIST tO 
STS  WhICh 
l  AS 7  INJUR 
XtMl  A/7 nt  P 
]Ut/ BLOOD 
1  r- S./P  40  ICC 
X/LT IOlOC-Y 
fcSSuKfc  RtS 


Ab7  1VJUR1FS  TO  7 HE  LUNGS  AFTER  BOMB 
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Abstract:  Literature  on  underwater  blast  effects  .in  man  and  animals 

was  reviewed  with  particular  reference  to  its  net  hoi  0.7  n 
r.i  a  t  h  o  n  h  u  s  i  o  1  o  a  u  a  n  d  t  h  e  r  a  p  u .  An  a  t  o  m  i  c  s  t  r  u  c  t  u  r  e  s  w  h  i  c  h  c:  o  n  t  a  i  n  a  i  r  ? 
i.e.?  lunosi  enteric  tract?  nasal  sinuses  and  middle  ear  were  found  to 
be  most  vulnerable  to  blast  iniuru.  An  historical  review  of 
therapeutic  procedures  used  in  the  treatment  of  blast  in  iu.ru  was  then 
presented .  Factors  found  to  be  of  area test  potential  benefit  in 
i  risrovino  the  dismal  survival  rate  of  underwater  blast  victims 
includes'  (1)  prevention  of  air  emboli?  (2)  maintenance  of  adeauato 
ventilation  and  respiration  and  <3>  time.'lu  suraical  repair  of  enteric 
tract  in  iuri.es . 
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Gians  and  symptoms?  Literature  surveys?  F'athcloau?  Cardiovascular 
system?  Hyperbaric  chambers?  Gas  embolism?  Gastroi ntestivial  system- 
Respiratory  system?  Mortality  rates.  J'liaanosis < Med.i cine )  ?  Anesthesia. 
First,  aid?  Gxuaen?  Blast  waves 
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Abstract,;  Underwater  blast  in  .mries?  at  increasi.no  r  a  acres  beyond  the 
lethal  2" one  from  small  charoes?  were  studied  usina  animals.  The  study 
was  conducted  in  an  artificial  pond  that,  measured  22''  iyy  ;  so  ft.  ot  its 
surface.  The  pond  was  30  ft  deep  over  .its  30-  bu  3  oo-f  t  center 

portion.  Sheep?  dorrs?  and  a  few  ronReus  were  exposed  to  the  blast 

oriented  vertically  in  the  water  (loner  axis  perpend.' tu  I  ar  to  the 
surface)  Most  were  exposed,  to  the  blast  at  i  ~ •+ +  dent hs •  t?  •ad-  ahov 
the  surface?  and  1.1.  mi.  ted  number  0+  .?—  and  10 ft  depths  F  ••'?.- lost  ve 

Ft  ?:.i  j”  jj .r. V' «-•  r> ?;% r-  I  j  1 1  ho  <"  Y't ,-j  i"* f  <  •’  •  •. :  ,  .  •  , 


•:<  lt».  Hi  1  chavoes  were  detonated  at  3  0- ft  depths;  The 

s  romersion-blast  in  iuri.es  were  of  minor  severity  and.  cons*  sied  ma.t  nl  •• 
of'  lunar  henorrhoaes  and  small  area#  of  contusions  in  r  he 
bistro: intestinal  t ra rt .  Tne  incidence  and  severity  of  the  in  mnes 
were  correlated  with  the  impulse  in  the  underwater  blast  wave  Based 
on  the  results  of  the  studio  a  safe  imtmlse  level  of  2  to  3  psa  nr-f" 
for  u  nor  o  tec  ted  swimmers.,  head  above  the  surface.-  was  nr  one  sea  Th  i.  =• 
safe  impulse  level  was  discussed  an  relation  to  the  underwater 
Mast-wave  parameters  in  the  test  pond  and  ex  is.;  tin  a  response  data  for 
p  e  r  s  o  n  vi  e  1  .  (Mod  ,i.  f  i  e  d  a  u  t.  h  o  r  a  b  s  t  r  a  <:  t ) 
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A b s t r a  c t  •  P r e s s u r e  p u  1.  s e  e x p e r i m e n t s  w e r e  conduc t e d  a n d  i n t  r a  c o r  p o r e a  3 
pressures  were  measured,  in  miniature  pigs  and.  albino  rats  using  a 
p r eviouslu  v a  1  i d a t e d  t e c h n i a u e  a n d  s p e c i a  1 1  u  s e  1  e c t e d  p r e s s u r e  r* r o b e s 
These  were  localized  for  the  miniature  pi as  in  the  esophagus.  the 
rectum,  and  in  the  musculature  of  the  back  and  the  thiah,  and  for  the 
albino  rats,  in  the  rectum.  The  clinical  symptoms  cause  by  pressure 
pulse  and  the  detected  morpholoaical  fa.nd.inas>  are  demonstrated  and 
verified  by  macro-  and  microscopical  photoaranhs.  The  variations  in 
the  results  of  these  experiments  are  discussed.  (Author') 

D  e  s  c  r  i  p  t  o  r  s  ■  *  P  r  e  s  s  u  r  e  a  r  a  d  i  e  n  t  s  •  *  R  a  t  s ,  *  S  w  i  n  e ,  D  i  a  e  s  t  i  v  e  s  u  s  t  e  m . 
Explosions.  Muscles,  Pressure  aaaes.  Pressure  measurements,  Rectum, 
Shock  waves 
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.  -  V  .* 

.-•mu*.*, utfe  .jv  -,2-mch  diapieter  shock  tube'.  Uhe' 'add  o  +  eacfi  ‘pdaV 'then 
was  given  intermittent  positive  pressure  re.~pi  ration  <  3  PF*R>  to*  2 
hours  with  100  percent  oxygen*  and  the  other  doer  was  maintained  on  :!.  oo 
percent  oxygen  for  4  hours  in  a  buperbaric  chamber  at.  o  chomber 
pressure  of  14  o .  s .  i  .  a .  *  after  which  she  was  aiven  If PR  with  ] 00 
tier  cent  oxuaen  for  2  hours.  The  mortal  itu*  tine  o+  death..  and 
inci device  of  arterial  air  embolism  in  these  two  groans  then  were 
compared  with  those  of  10  untreated  control  am. no. is  that  previous!  u 
had  been  exposed  to  airblast  in  the  same  wan  as  those  in  the  treatment 
groups.  The  mortal itu  was  60  percent  i  vt  t  he  untreated  control  a rou y.:>* 
BO  percent  in  the  immediate  IF'F’R  group .  and  50  percent  in  the  del  oued 
IPPR  treatment  group.  There  was  one  case  of  air  embolism  ( 1 4-rs.iviuto 
fatality)  in  the  untreated  control  group*  three  cases  of  air  embolism 
in  the  immediate  IPPR  group*  and  none  in  the  del  a  tied  IPPR  grout*.  The 
mean  survival  time  for  the  fatalities  was  12.4  hours  for  the  untreated 
control  group*  2.3  hours  for  the  immediate  IPPR  group*  and.  9  9  hours 
for  the  del  aired  IPF'R  group.  Thus*  the  results  indicate  that  the  use 
of  IF'F'R  immediate lir  following  blast  in  luru  mat./  result  i  n  an  increase 
in  the  incidence  of  air  embolism*  increase  in  mart a lit  to  and  a 
reduction  in  survival  time-  whereas*  when  used  after  a  del  air  of  4 
hours*  IF'F'R  resulted  in  neither  an  increase  in  incidence  of  air 
embolism  nor 'in  mortal itu  but  did  result  in  a  shortening  of  survival 
time.  (Author) 

Descriptors :  (* Pressure  breathing,  Theraou)  *  (i'Gas  embolism*  Pressure 
breathing)*  (fcBlast*  Gas  embolism)*  OKLungs*  Explosion  effects  >  * 
Wounds  +  In  luries*  Dxirgen*  Mortality  rates*  Survival*  Shock (  Pathology) 
*  Aviation  medicine*  Do as*  Experimental  data 
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Abstract:  Guinea  pigs  avid  rabbits  were  exposed  to  lethal  reflected 
pressures  in  an  air-driven  shock  tube  avid  were  subseauentlu  treated  in 
a  hyperbaric  chamber  in  which  the  oxugen  tension  (P02)  and  chamber 
pressure  were  .independent lu  varied.  Treatments  involving  increases  in 
P02  resulted  in  increased  survival  times  of  guinea  Digs  whereas 
pres  sari  rat. ion  for  30  minutes  at  36  or  72  p.s.i.g.  with  the  F'02 
retained  at  the  normal  ambient  level  bu  use  of  an  N2-a.tr  mixture  had 
no  detectable  effect  on  survival  times  of  the  avtintals.  To  study  the 
effects  of  prolonged  hyperbaric  oxygenation  in  treatment-  of  blast 
ivtiuru*  guinea  pigs  and  rabbits  were  treated  on  a  2 9 -hour  schedule 
having  art  initial  3-hour  hold-time  at  the  pressure-treatment  level 
followed  bu  26  hours  for  decompression  In  rabbits*  an  .initial  PQ2  of 
.1.7.5  p.s.i .a  *  achieved  either  bu  air  pressure  at  72  p.s.i.g.  or  bu 
* .  ressut-i? a t i o n  to  1.5  p  s.i.g.  wi Jh  65— per-cent  0? «  35— percent  *; ? * 
resulted  in  tall  survival  and  recovery  o->"  all  treated  animals.  In 
guinea  pigs*  treatment  with  100-T<ev-cent  02  at  5.  5  t*  «:  .  i  a  (PG2  -  1.7  ■* 
r«.  s.  i.  a.  )  or  at  12  p.s.i.  a  (PG?  -  2y  p.s.i  a.  >  resulted  in  inrrease*.4 
r  Mi'  viva  I  times  with  no  .increase  in  overall  survival  and  recovery  tv* 
the  firs*  case  and  significant  lu  increased  survival  and  recover** 
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Abstract :  The  blast  effects  in  rectangular  two-nan  foxholes  were 
evaluated  using  sheen.  There  were  two  open  foxholes  at  around  ranges 
of  560,  650,  830,  940,  and  1,300  feet  from  a  500-ton  TNT  char  ore. 
Decause  of  an  anomalous  detonation,  pressures  measured  ad  lacent  to  the 
•foxhole  lauout  were  significant, lu  below  those  predicted  Moreover, 
luminous  lets  emanatina  from  the  fireball  produced  shock  waves  that 
preceded  the  main  shock.  This  aave  rise  to  a  blast  wave  with  double 
shocks  known  oeueral  lu.  to  be  less  damaaina  to  bio  load,  cal  s  us  terns.  All 
the  sheep  survived  the  blast.  At  the  560-  and  650 -foot  rartaes  <37  and 
21  o.s.i  )  some  of  the  sheen  sustained  sliaht  amounts  of  pul  mortar  u 
hemorrhaae.  In  addition,  th.eu  exhibited  a  hi  ah  incidence  of  eardrum 
r u p t u r e  of  a  s e v e r e  f orm.  < A u t h o r ) 
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■.in  mrU)  the  pathol  ocru,  and  cunicai  consideration*  are  discuss*? 
discussion  and  criticise  is  presented  of  the  various  formula 
da maae  rancre.  Much  of  the  material  is  supported  with  reference 
b  a  t  h  a  n  i  m  a  1  a  n  d  h  u  m  a  n  d  a  t,  a  .  (A  u  t  h  o  r  ) 
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Abstract:  A  us  den  vorlieerenden  Untersuchunasergebnissen  geht  hervor, 
d  a  s  s  bei  Z  w  e  r  or  s  c  h  w  e  i.  n  e  n  n  a  c  h  S  c  h  a  d  i.  a  u  n  a  d  u  r  c  h  D  r  u  c  k  w  e  1 1  e  n  s  t  o  s  s  d  .1  e 
Uberdruckbehandluna  in  Verbinduna  nit  diuretischen  Hassnahnen*  wie  mit 
der  Verabriechung  von  , Lasix' (Fursemide),  ohne  Zweifel  a  Is  eine 
m o a  1  i  ch e  ■  "T h  e r a p .i e  F  o r m  i h r e  B e d e u t u n c  h.  a t .  Be r  k r i t .i s  c h e  e r a  1  e i  c h  v o n 
Thera pie-  und  Kontr  >11 tier en  in  Gruppen,  die  21 — 35  Tags  und  36 — 60 
f age  uherlebt  ha ben,  hat  eraeben,  dass  bei  den  behandelten 
Zwergschweinen  nach  an  sich  schneller  Resorption  der  Blutunoen,  die 
reparativen  Vorgange  in  den  Lunaeri  spa  ter  einsetzen  und  weriiger 
ausgepracrt  sind.  Die  Behan dlung  hat  bei  den  druckst ossa eschadig ten 
Versuchstieren  die  moroholooisch  nachaeuiesenen  Spa tver under ungen  in 
den  Lungen  als  Folae  der  Blutunaen  nicht  verhindern,  aber  doch 
zum.i  ndest  einsch  ranker  konnen.  Da  bei  muss  man  der  Ga.be  von  ,  Lasix'  mit 
der  prompter,  d iuretischen  Uirkuna  einen  besonderen  theraoeutischen 
Fffe k t  b e i m e s s e n .  (A u t h o r ) 
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Abstract :  The  pattern  of  recover;.-  e  f  The  respi  ratoi'u  sustem  f  ron  blast 
i  r.  luru  was  invest. icrated  in  sheen  exposed  to  over pro  s  sur  e  i  >•  a  short 
tube  .  Measurements  of  the  pH  and  biooa  aa«  t  ensions  del e-v-mriut  .1  on- 

of  the  venous-admixture  (Os/!?/  and  the  al.veoiarav-i.ey  ini  o  •  inr.-n 
Gradient-  <A-a)C2  were  conducted  be  fore  and  at  .interval  s  ut>  to  J  32  dans 
t  olloi/.una  in  luv'i.r  There  was  an  immediate  narked  increase  i-»  U-  1  , 
reduction  in  Po02,  and  a  moderate  increase  in  (A- a  >02?  with  v:  ru 
1  ittl*  chanae  in  the  pH  or  PC02  of  +he  arterial  Mood. .  ( (>>■-  great  e- f 

r  ecoveru  was  evident  within  24  hours  with  further  aroduo  .  ii-.prov- -m- -  - 
seen  2.  7 ,  14/  and  21  da  us  after  exposure.  After  the  21st  dan.  wort  o'" 

the  animals  exhibited  virtual  complete  recovery  of  the  functional 
of f iciencu  of  the  pulmonaru  sustem  as  tested  at  rest.  <  Author' 
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Abstract ••  Shock  tubes  and  hi. ah  explosives  were  used  to  Drodu.ce  Mast 
waves  of  various  pressure-time  cat terns  in  order  to  studu  their 
b  i  o  1  o  a  i  c  a  1  e  f  f  &  c  t  s .  D  a  t  a  o  b  t  a  i  n  e  d  f  r  o  m  t  h  e  s  e  e  x  p  e  r  i  n  e  n  t  s  s  h  o  w  e  d  t  h  a  t  - 
aaainst  a  ref lectina  surface,  the  LD50  reflected,  pressure  for  am: 
aiven  species  remained  fairlu  constant  at  the  'ionaer '  durations  and 
then  rose  share lu  at  the  "shorter"  times.  For  do as  and  aoats,  'lone 
durations  were  beuond  20  msec  and  for  mice,  rats,  eruiviea  Digs,  and. 
rabbits,  beuond  1  to  3  msec  At  the  'shorter '  durations,  response 

depended  to  a  erroat  extent  on  the  impulse-  and  on  peak  pressure  +ov 
the  "longer"  pulses.  Hi  a  her  reflected  pressures  can  be  withstood  if 
animals  are  located  beuond  a  certain  distance  from  the  re* lectina 
surface  where  then  receive  the  incident  and  reflected  pressures  in  two 
steps,  separated  bu  a  aiven  time-interval.  In  frees  tv*  earn  exposures  to 
air  blast-  orientation  was  sianif icant .  Animals  suspended  vertical lu 

or  prone . side -on  showed  a  lower  to  lev*  a  nee  to  blast  waves  of  a  aiven 

intensity  or  at  a  aiven  ranae  than  those  end-on  because  the  dunam?,  r 
pressure  appeared  to  add.  to  their*  side-on  pressure  dose.  Except  *  o~ 
oardrur-i  rupture  and  sinus  hemorrhaae,  animals  exhibited  o  rewari  a  M  a 
tolerance  to  "slow'-risina  blast  pressures  without  the-  presence 
shod:  fronts.  The  lunar s  are  considered  the  critical  target  ov*gans  in 
b  1  a  s  t  e  f  f  e  c  t  s  s  t  u  d.  i  e  s .  ( A  '.j  t  h  a  r ) 
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Abstract:  Blast  iniuru  is  a  Comdex  and  veru  hazardous  phenomenon  to 

the  biolocric  target.  Tocrether  with  effects  of  thermal  radiations 
■f rr  i  modern  nuclear  weapons.*  blast  in  luru  (direct  and  indirect  ) 
an  nears  to  be  accountable  for  the  vast  bull  of  earlu  deaths  avid 
casualties  .in  nuclear  explosions.  This  article  has  attempted  to 
summarize  the  important  clinical/  oh  usio  logic.-  and  pathologic 
information  concerning  the  effects  o fdjiie 1  a s t  in  mru  on  the 
biologic  subiect.  Certain  feat  u  u'(sr  n  rv  a  s  i.zed  in  order  to 

assist  the  clinical  medical  officer  towards  proper  management  of 
casualties.  A  brief  description  of  pulmonary  sequelae  of  blast  i niu.ru 
is  included  for  completeness.  (Author) 
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Abstract:  The  purpose  of  this  study  was  to  predict  the  probability  of 

impact  injuries  due  to  whole-body  translation  by  airblast  as  a 
function  of  yield  and  ground  range.  Predictions  were  made  for 
personnel  in  different  orientations  in  open  terrain  and  near 
structural  complexes.  A  mathematical  model  was  used  to  calculate  the 
time-displacement  history  of  personnel  from  considerations  of 
aerodynamic  drag  and  ground  friction.  Predicted  values  of  maximum 
velocity*  displacement  at  maximum  velocity*  and  total  displacement 
were  tabulated  for  1224  exposure  conditions.  Biological  criteria  were 
presented  which  indicated  that  personnel  subjected  to  decelerative 
fumbling  over  open  terrain  can  tolerate  much  higher  velocities  than 
personnel  impacting  a  nonuielding*  flat  surface  at  normal  incidence. 
Methods  for  extending  the  presented  results  to  other  exposure 
conditions  were  discussed. 


exposure 


Descriptors:  ^Impact  shock*  *Airburst*  *Nuclear  warfare*  ^Tumbling* 
founds  and  in  luries,  Yield (Nuclear  explosions)*  Range (Distance ) > 
Mathematical  models*  Casualties*  Aerodynamic  drag*  Velocity*  Blast 
caves*  Displacement*  Humans 

Identifiers:  NTISDODXA*  NTISDODSD 

AD— AOS?  785/ OST  NTIS  Prices:  PC  A03/MF  A01 
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